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SECTION  1.  GENERAL  INFORMATION 


1.1  BACKGROUND 

Technologies  under  development  for  the  detection  and  discrimination  of  unexploded 
ordnance  (UXO)  require  testing  so  that  their  performance  can  be  characterized.  To  that  end, 
Standardized  Test  Sites  have  been  developed  at  Aberdeen  Proving  Ground  (APG),  Maryland, 
and  U.S.  Army  Yuma  Proving  Ground  (YPG),  Arizona.  These  test  sites  provide  a  diversity  of 
geology,  climate,  terrain,  and  weather  as  well  as  diversity  in  ordnance  and  clutter.  Testing  at 
these  sites  is  independently  administered  and  analyzed  by  the  Government  for  the  purposes  of 
characterizing  technologies,  tracking  performance  with  system  development,  comparing 
performance  of  different  systems,  and  comparing  performance  in  different  environments  (app  E, 
ref  1). 


The  Standardized  UXO  Technology  Demonstration  Site  Program  is  a  multiagency 
program  spearheaded  and  funded  by  the  Environmental  Securities  Technology  Certification 
Program  (ESTCP),  the  Strategic  Environmental  Research  and  Development  Program  (SERDP). 
The  U.S.  Army  Aberdeen  Test  Center  (ATC)  provides  programmatic  and  field  support  for 
technology  demonstration  and  evaluation,  and  maintains  a  repository  of  inert  munition  items 
available  to  the  UXO  community.  The  U.S.  Army  Environmental  Command  maintains  the 
Standardized  UXO  Technology  Demonstration  Site  Program  web  page 
(http://aec.army.mil/usaec/technology/uxoOEhtml),  which  contains  program  information,  vendor 
demonstration  instructions  and  copies  of  all  published  vendor  demonstration  scoring  records. 

1.2  SCORING  OBJECTIVES 

The  objective  in  the  Standardized  UXO  Technology  Demonstration  Site  Program  is  to 
evaluate  the  detection  and  discrimination  capabilities  of  a  given  technology  under  various  field 
and  soil  conditions.  Inert  munitions  and  clutter  items  are  positioned  in  various  orientations  and 
depths  in  the  ground. 

The  evaluation  objectives  are  as  follows: 

a.  To  determine  detection  and  discrimination  effectiveness  under  realistic  scenarios  with 
various  targets,  geology,  clutter,  density,  topography,  and  vegetation. 

b.  To  determine  cost,  time,  and  workforce  requirements  to  operate  the  technology. 

c.  To  determine  the  demonstrator’s  ability  to  analyze  survey  data  in  a  timely  manner  and 
provide  prioritized  Target  Eists  with  associated  confidence  levels. 

d.  To  provide  independent  site  management  to  enable  the  collection  of  high  quality, 
ground-truth  (GT),  geo-referenced  data  for  post-demonstration  analysis. 
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1,2.1  Scoring  Methodology 

a.  The  scoring  of  the  demonstrator’s  performance  is  conducted  in  two  stages:  response 
stage  and  discrimination  stage.  For  both  stages,  the  probability  of  detection  (Pa)  and  the  false 
alarms  are  reported  as  receiver-operating  characteristic  (ROC)  curves.  False  alarms  are  divided 
into  those  anomalies  that  correspond  to  emplaced  clutter  items,  measuring  the  probability  of 
clutter  detection  (Pcd)  or  the  probability  of  false  positive  (Pfp).  Those  that  do  not  correspond  to 
any  known  item  are  termed  background  alarms.  The  background  alarms  are  addressed  as  either 
probability  of  background  alarm  (Pba)  or  background  alarm  rate  (BAR). 

b.  The  response  stage  scoring  evaluates  the  ability  of  the  system  to  detect  emplaced 
targets  without  regard  to  ability  to  discriminate  munitions  from  other  anomaly  sources.  For  the 
blind  grid  response  stage,  the  demonstrator  provides  a  target  response  from  each  and  every  grid 
square  along  with  a  threshold  below  which  target  responses  are  deemed  insufficient  to  warrant 
further  investigation.  This  list  is  generated  with  minimal  processing  and,  since  a  value  is 
provided  for  every  grid  square,  includes  amplitudes  both  above  and  below  the  system  noise  level. 
For  the  open  field,  the  demonstrator  provides  a  list  of  all  anomalies  deemed  to  exceed  a 
demonstrator  selected  target  detection  threshold.  An  item  (either  munition  or  clutter)  is  counted 
as  detected  if  a  demonstrator  indicates  an  anomaly  within  a  specified  distance  (Halo  Radius 
(Rhaio))  of  a  ground  truth  item. 

c.  The  discrimination  stage  evaluates  the  demonstrator’s  ability  to  correctly  identify 
munitions  as  such  and  to  reject  clutter.  For  the  blind  grid  discrimination  stage,  the  demonstrator 
provides  the  output  of  the  discrimination  stage  processing  for  each  grid  square.  For  the  open 
field,  the  demonstrator  provides  the  output  of  the  discrimination  stage  processing  for  anomaly 
reported  in  the  response  stage.  The  values  in  these  lists  are  prioritized  based  on  the 
demonstrator’s  determination  that  a  location  is  likely  to  contain  munitions.  Thus,  higher  output 
values  are  indicative  of  higher  confidence  that  a  munitions  item  is  present  at  the  specified 
location.  For  digital  signal  processing,  priority  ranking  is  based  on  algorithm  output.  For  other 
discrimination  approaches,  priority  ranking  may  be  based  on  rule  sets  or  human  judgment.  The 
demonstrator  also  specifies  the  threshold  in  the  prioritized  ranking  that  provides  optimum 
performance,  (i.e.,  that  is  expected  to  retain  all  detected  munitions  and  reject  the  maximum 
amount  of  clutter). 

d.  The  demonstrator  is  also  scored  on  efficiency  and  rejection  ratios,  which  measure  the 
effectiveness  of  the  discrimination  stage  processing.  The  goal  of  discrimination  is  to  retain  the 
greatest  number  of  munitions  detections  from  the  anomaly  list,  while  rejecting  the  maximum 
number  of  anomalies  arising  from  nonmunitions  items.  Efficiency  measures  the  fraction  of 
detected  munitions  retained  after  discrimination,  while  the  rejection  ratio  measures  the  fraction 
of  false  alarms  rejected.  Both  measures  are  defined  relative  to  the  maximum  number  of 
munitions  detectable  by  the  sensor  and  its  accompanying  clutter  detection/false  positive  rate  or 
BAR. 
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e.  Based  on  configuration  of  the  GT  at  the  standardized  sites  and  the  defined  scoring 
methodology,  in  some  cases,  there  exists  the  possibility  of  having  anomalies  within  overlapping 
halos  and/or  multiple  anomalies  within  halos.  In  these  cases,  the  following  scoring  logic  is 
implemented; 

(1)  In  situations  where  multiple  anomalies  exist  within  a  single  Rhaio,  the  anomaly  with 
the  strongest  response  or  highest  ranking  will  be  assigned  to  that  particular  GT  item.  If  the 
responses  or  rankings  are  equal,  then  the  anomaly  closest  to  the  GT  item  will  be  assigned  to  the 
GT  item.  Remaining  anomalies  are  retained  and  scored  until  all  matching  is  complete. 

(2)  Anomalies  located  within  any  Rhaio  that  do  not  get  associated  with  a  particular  GT 
item  are  excess  alarms  and  will  be  disregarded. 

f  In  some  cases,  groups  of  closely  spaced  munitions  have  overlapping  halos.  The 
following  scoring  logic  is  implemented  (App  A,  fig.  A-1  through  A-9): 

(1)  Overall  site  scores  (i.e.,  Pd)  will  consider  only  isolated  munitions  and  clutter  items. 

(2)  GT  items  that  have  overlapping  halos  (both  munitions  and  clutter)  will  form  a  group 
and  groups  may  form  chains. 

(3)  Groups  will  have  a  complex  halos  composed  of  the  composite  halos  of  all  its  GT 

items. 

(4)  Groups  will  have  three  scoring  factors:  groups  found,  groups  identified,  and  group 
coverage.  Scores  will  be  based  on  1 ;  1  matches  of  anomalies  and  GT. 

(a)  Groups  Found  (Found):  the  number  of  groups  that  have  one  or  more  GT  items 
matched  divided  by  the  total  number  of  groups.  Demonstrators  will  be  credited  with  detecting  a 
group  if  any  item  within  the  group  is  matched  to  an  anomaly  in  their  lists. 

(b)  Groups  Identified  (ID):  the  number  of  groups  that  have  two  or  more  GT  items 
matched  divided  by  the  total  number  of  groups.  Demonstrators  will  be  credited  with  identifying 
that  a  group  is  present  if  multiple  items  within  the  composite  halo  are  matched  to  anomalies  in 
their  lists. 

(c)  Group  Coverage  (Coverage):  the  number  of  GT  items  matched  within  groups  divided 
by  the  total  number  of  GT  items  within  groups.  This  metric  measures  the  demonstrator  accuracy 
in  determining  the  number  of  anomalies  within  a  group.  If  five  items  are  present  and  only  two 
anomalies  are  matched,  the  demonstrator  will  score  0.4.  If  all  five  are  matched,  the  demonstrator 
will  score  1.0. 

(5)  Location  error  will  not  be  reported  for  groups. 
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(6)  Demonstrators  will  not  be  asked  to  eall  out  groups  in  their  seoring  submissions.  If 
multiple  anomalies  are  indieated  in  a  small  area,  the  demonstrator  will  report  all  individual 
anomalies. 

(7)  Exeess  alarms  within  a  halo  will  be  disregarded. 

g.  All  seoring  faetors  are  generated  utilizing  the  Standardized  UXO  Probability  and  Plot 
Program,  version  4. 

1,2,2  Scoring  Factors 

Factors  to  be  measured  and  evaluated  as  part  of  this  demonstration  include: 

a.  Response  stage  ROC  curves: 

(1)  Probability  of  detection 

(2)  Probability  of  clutter  detection  (Pcd)- 

(3)  Background  alarm  rate  (BAR"^®*)  or  probability  of  background  alarm  (Pba'^*^'^). 

b.  Discrimination  stage  ROC  curves: 

(1)  Probability  of  detection 

(2)  Probability  of  false  positive  (Pfp). 

(3)  Background  alarm  rate  (BAR‘*'®®)  or  probability  of  background  alarm  (Pba^'*^^)- 

c.  Metrics: 

(1)  Efficiency  (E). 

(2)  False  positive  rejection  rate  (Rfp). 

(3)  Background  alarm  rejection  rate  (Rba). 

d.  Other: 

(1)  Probability  of  detection  by  size,  depth,  and  density. 

(2)  Classification  by  type  (i.e.,  20-,  40-,  105-mm,  etc.). 

(3)  Eocation  accuracy  for  single  munitions. 
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(4)  Equipment  setup,  calibration  time,  and  corresponding  worker-hour  requirements. 

(5)  Survey  time  and  corresponding  worker-hour  requirements. 

(6)  Reacquisition/resurvey  time  and  worker-hour  requirements  (if  any). 

(7)  Downtime  due  to  system  malfunctions  and  maintenance  requirements. 
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SECTION  2.  DEMONSTRATION 


2.1  DEMONSTRATOR  INFORMATION 

2.1.1  Demonstrator  Point  of  Contact  IPOC)  and  Address 

POC:  Ted  Asch 

(303)  236-2489 

Address:  U.S.  Geologieal  Survey 
Denver  Federal  Center 
Bldg.  20,  MS-964 
Denver,  CO  80225-0046 

2.1.2  System  Description  (provided  by  demonstrator) 

The  ALLTEM  is  an  ‘on-time’  time-domain  eleetromagnetie  induetion  system  that  eonsists 
of  exeiting  and  deteeting  3-oomponent  fields  using  multiple  Tx  and  Rx  eoils.  The  triangle  eurrent 
exeitation  waveform  (pulse  rate  90Hz)  provides  immediate  visual  separation  between  ferrous  and 
non-ferrous  metal  objects.  The  ALLTEM  records  data  to  late  times  which  helps  suppress  the 
geologic  response  relative  to  the  UXO  response.  The  system  is  pulled  by  a  small  Kubota  tractor 
with  a  small  2  kW  generator  at  the  front  (see  photo  below).  The  ALLTEM  sensor  is  a  1 -meter 
cube  that  sits  in  a  cart  and  has  a  minimum  height  above  the  ground  of  about  6  inches  that  can  be 
raised  up  an  additional  6  inches  to  traverse  over  surface  obstacles.  Two  LEICA  GPS  1200  RTK 
rover  systems  provide  the  sensor  location  and  tractor  position  into  a  USGS-developed  survey 
navigation  program.  Survey  traverses  will  have  0.5  meter  separation  with  a  data  density  of 
approximately  15  to  20  cm  (traveling  at  a  nominal  speed  of  1.0  m/sec  with  a  sampling  cycle  rate 
for  each  Tx  coil  of  approximately  300  ms). 
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Figure  1.  Demonstrator’s  system,  EM61  MKII/towed  array. 


2,1,3  Data  Processing  Description  (provided  by  demonstrator) 

Target  selection  criteria:  This  section  will  detail  the  target  selection  criteria  and  the  data 
required  to  implement  the  criteria  by  answering  the  following  questions: 

a.  What  kind  of  pre-processing  (if  any)  is  applied  to  the  raw  data  (e.g.,  filtering,  etc.)? 
ALLTEM  preprocessing  is  a  batch  process  of  all  binary  waveform  survey  data  via  a  EABVIEW 
program  that  performs  background  subtraction,  low-pass  and  band-width  filtering,  determination 
of  ferrous/nonferrous/mixed  composition,  and  then  exports  an  ASCII  file  containing  data  at 
16  time  gates  along  the  waveform. 

b.  What  is  the  format  of  the  data  both  pre  and  post  processing  of  the  raw  data  (e.g., 
ASCII,  binary,  etc.)?  The  original  EABVIEW  acquisition  data  consists  of  binary  waveform  files 
with  ASCII  headers.  There  is  one  file  per  configuration.  These  are  converted  in  the  EABVIEW 
preprocessing  program  to  an  ASCII  format  that  is  carried  throughout  the  rest  of  the  processing 
and  analysis. 

c.  What  algorithm  is  used  for  detection  (e.g.,  peaks  of  signal  surpassing  threshold,  etc.)? 
In  2008  we  have  migrated  all  our  processing  and  analysis  software  to  work  within  the  Geosoft 
Oasis  Montaj  platform.  Once  the  data  is  imported  into  Oasis,  an  area  that  is  deemed  to  be  target 
free  is  designated.  This  area  forms  the  threshold  basis  on  which  a  statistical  analysis  is  performed 
using  the  “R  Project  for  Statistical  Computing”  (http://www.r-project.org)  statistics  package, 
version  2.10.1.  Wilkes-Shapiro  and  T-tests  characterize  the  acquired  data  and  then  Blakely 
peakedness  tests  are  performed  to  designate  the  locations  of  the  potential  targets.  This  is  all  done 
automatically  for  all  19  AEETEM  receiver  configurations. 
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d.  Why  is  this  algorithm  used  and  not  others?  We  use  the  ealeulated  statisties  for  both 
pieking  out  targets  and  as  part  of  the  elassifieation  analysis  at  the  end  of  the  proeess. 

e.  On  what  principles  is  the  algorithm  based  (e.g.,  statistical  models,  heuristic  rules, 
etc.)?  As  just  mentioned,  the  algorithm  is  based  on  a  statistieal  analysis  of  the  aequired  data. 

f  What  tunable  parameters  (if  any)  are  used  in  the  detection  process  (e.g.,  threshold  on 
signal  amplitude,  window  length,  filter  coefficients,  etc.)?  Tunable  parameters  inelude  the 
baekground  threshold  level,  the  number  of  standard  deviations  away  from  the  target  threshold 
used  to  determine  signal  levels,  the  seareh  radius  around  eaeh  seleeted  target  (used  for  merging 
multiple  targets  at  same  loeation  from  different  reeeiver  polarizations),  the  areas  of  what  are 
eonsidered  to  be  statistieally  ‘signifieant’  data  for  a  partieular  target,  and  analytie  signal 
ealeulations  for  eertain  reeeiver  polarizations. 

g.  What  are  the  final  values  of  all  tunable  parameters  for  the  detection  algorithm?  The 
final  values  for  the  tunable  parameters  will  be  determined  by  the  data  in  the  field.  The 
baekground  threshold  values  will  be  determined  independently  for  eaeh  area  surveyed.  The 
seareh  radius  will  be  determined  by  the  largest  target  deteeted  in  eaeh  survey  area. 

Parameter  estimation:  This  seetion  should  inelude  the  details  of  whieh  parameters  will 
be  extraeted  from  the  sensor  data  for  eaeh  deteeted  item  for  eharaeterization.  Please  answer  the 
following  questions: 

a.  Which  characteristics  will  be  extracted  from  each  detected  item  and  input  to  the 
discrimination  algorithm  (e.g.,  depth,  size,  polarizability  coefficients,  fit  quality,  etc.)? 
Charaeteristios  extracted  for  each  detected  item  include  inferred  composition 
(ferrous/nonferrous/mixed),  horizontal  location  and  depth,  azimuth,  inclination,  magnetic 
polarizability  coefficients,  and  the  ratio  of  polarizability  coefficients. 

b.  Why  have  these  characteristics  been  chosen  and  not  others  (e.g.,  empirical  evidence  of 
their  ability  to  help  discriminate,  inclusion  in  a  theoretical  tradition,  etc.)?  We  have  determined 
empirically  from  previous  surveys  and  by  models  that  these  characteristics  (composition, 
polarizabilities,  ratios  of  polarizabilities)  have  proven  effective  at  discriminating  UXO  versus 
clutter  versus  blank  holes. 

c.  How  are  these  characteristics  estimated  (e.g.,  least-mean-squares  fit  to  a  dipole  model, 
etc.),  include  the  equations  that  are  used  for  parameter  estimation? 

The  ALLTEM  response  is  modeled  using  a  conductive,  magnetic,  and  optionally,  a  viscous 
magnetic  earth.  The  integral  expressions  for  the  magnetic  fields  generated  by  infinitely  small 
vertical  magnetic  dipoles  (VMDs)  and  horizontal  magnetic  dipoles  (HMDs)  over  a  conductive 
magnetic  half  space  are  given  by  Ward  and  Hohmann  (1988).  These  integrals  are  numerically 
evaluated  using  Anderson’s  (1979)  fast  Hankel  transform  to  produce  Green’s  functions  for  the 
magnetic  field  due  to  a  dipole  over  a  half-space.  The  ground  response  at  a  given  receiver  coil 
due  to  excitation  by  a  given  transmitter  coil  is  calculated  by  numerically  evaluating  the  following 
integrals. 
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(1) 


V^{o})  =  -i(oT{(o)  j  ^hi^-G  {(T,/u,r,r')drdr 

RXloop  TXloop 


where  T(co)  is  the  Fourier  transform  of  the  transmitter  exeitation  (a  triangle  wave), 
G  (cr,//,r,r')is  the  dyadie  Green’s  funetion  for  the  fields  above  the  surfaee  due  to  a  dipole 
moment  above  the  surfaee,  a  is  the  ground  eonduetivity,  ji  is  magnetie  permittivity  (possibly 
eomplex  and  frequeney  dependent),  co  is  radian  frequeney,  Vrx  is  the  voltage  at  the  reeeiver  eoil, 
and  are  the  normals  to  the  planes  eontaining  the  transmitting  and  reeeiving  eoils,  and  the 
integrals  are  over  the  area  of  the  transmitting  and  reeeiving  eoils.  This  proeedure  is  used  to 
ealeulate  the  voltage  at  eaeh  reeeiving  eoil  in  the  gradiometer,  and  the  simulated  result  is  the 
differenee  between  these  eoil  voltages.  Finally,  the  frequency-domain  coil  voltage  is  converted 
into  a  time-domain  signal  using  a  numerical  Fourier  transform  (FFT). 

Modeling  the  response  of  a  conductive  permeable  spheroid  is  slightly  more  involved. 
First,  the  magnetic  field  below  the  surface  is  calculated  at  the  center  of  the  target  using  (2).  Here 
the  dyadic  Green’s  function  is  now  for  the  fields  below  the  surface  due  to  a  dipole  moment 
above  the  surface. 

//o(r,®)  =  r(®)  ^h,^-G^{(j,p,r,r)dr  (2) 

TXloop 


Next,  the  equivalent  induced  dipole  moment  m(t)  for  the  target  is  calculated, 

m{Q))  =  R^  ■M{Q))-R-Ho{r,Q))  (3) 

where  M{q))  is  the  diagonal  frequency-dependent  magnetic-polarizability  tensor.  Calculations 
for  specifying  this  tensor  for  conductive  permeable  prolate  and  oblate  spheroids  are  given  in 
Smith  and  Morrison  (2006).  is  a  rotation  matrix  that  converts  the  magnetic  field  to  target 
spheroid  centric  coordinates.  Note  that  for  a  given  incident  field,  only  a  single  component  of  the 
polarizability  tensor  is  excited.  A  fully  polarimetric  instrument  switches  through  different  coil 
configurations  resulting  in  incident  fields  with  components  in  all  possible  directions.  The  result 
is  that  all  components  of  the  polarizability  tensor  are  excited,  which  provides  more  information 
about  the  target. 

The  final  step  is  to  calculate  the  fields  at  the  receiver  coils  due  to  the  induced  dipole 
moments  at  the  target.  Using  the  reciprocity  theorem,  we  find  that  the  voltage  induced  in  the 

receiver  coil  F„  due  to  the  induced  target  dipole  moment  in  is  related  to  magnetic  field  at 
the  target  due  to  a  current  in  the  receiver  coil 


V^^{co)  =  -iQ} 


I 


rx 


(4) 
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Equation  2  is  used  to  calculate  the  magnetic  field  at  the  target,  only  now  the  integral  is  over  the 
area  of  the  receiver  coil.  As  before,  this  procedure  is  used  to  calculate  the  voltage  at  each 
receiving  coil  in  the  gradiometer,  and  the  simulated  result  is  the  difference  between  these  coil 
voltages.  Finally  the  frequency-domain  coil  voltage  is  converted  into  a  time-domain  signal  using 
a  numerical  Fourier  transform  (FFT). 

Using  Equations  1-4  and  Smith  and  Morrison’s  (2006)  spheroid  response,  the  AEFTEM 
response  to  the  earth  and  the  target  are  modeled  separately,  and  then  summed  for  the  total 
response.  Interactions  between  the  target  and  surrounding  medium  are  neglected  (i.e..  Bom 
approximation).  The  magneto-static  response  due  to  magnetization  of  a  permeable  target  is 
manifest  as  a  square  wave  response,  with  the  decaying  electro -dynamic  response  due  to  induced 
eddy  currents  superimposed.  The  square  wave  magneto-static  response  is  absent  for  the 
non-permeable  target.  Note  also  that  the  electro-dynamic  eddy  current  decay  lasts  longer  for  the 
permeable  target. 

When  processing  large  amounts  of  data,  it  is  desirable  to  have  a  fast  inversion  algorithm, 
which  in  turn  requires  a  fast  forward  model.  Calculating  the  Green’s  function  in  equations  1  and 
4  is  computationally  expensive,  and  is  not  needed  is  most  cases  (i.e.,  when  viscous  magnetic  soil 
is  not  present).  A  faster  approach  to  calculating  the  magnetic  fields  is  to  use  the  static 
Biot-Savart  Faw  for  free  space.  This  works  especially  well  for  the  AEFTEM  system  since  all 
coils  have  a  square  shape.  The  magnetic  field  at  r  is  the  sum  of  the  fields  produced  by  the  four 

straight  wire  segments,  each  with  current/^ ,  length  2L,  and  centered  atrj  (the  midpoint  of  the 
wire  segment): 


(^^2  ^  ~  -^n  )  (^1  ^  ~  -^n  ) 


(V) 


2 

2  1  , 

2 

L^  +  2L 

Z 

+ 

P 

+ 

z 

,  d,=JL^-2L 

Z 

+ 

P 

+ 

z 

^n=K-K^  Pn=K-Z„. 


(8-9) 

(10-12) 
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Smith  et  al.  (2004)  presented  a  simple  parametrie  form  for  estimating  the  time-domain 
B  field  response  of  a  eonduetive  permeable  sphere  due  to  a  step  funetion  exeitation.  This  is  also 
the  eleetro-dynamie  ALLTEM  response  sinee  its  exeitation  is  integral  of  the  step  (triangle  wave), 
it  uses  dB/dt  reeeivers,  and  the  integral  and  derivative  operations  eaneel  eaeh  other.  The  form  is  a 
good  approximation  to  the  early-time,  intermediate-time,  and  late-time  portions  of  the  response. 
The  simple  models  for  permeable  and  non-permeable  spheres  reduee  to: 


M,,Uereit,K) 


9/u^ 

3  2(//,+2) 


1  + 

V 


(13) 


a=1.38r„  ^ 

{Sl ,  =  ;z- +  arotan((^j//^  -^i)/(A.  + 

=  To ,  =  07/,//oi?7 ((a.  +  2)(a.  “  1))  • 


(14-16) 

(17-18) 

(19-20) 


We  eombined  this  form  with  Smith  and  Morrison’s  (2006)  approximation  of  the  polarizability 
tensor  for  a  prolate  spheroid  in  terms  of  that  of  a  sphere.  For  the  element  along  the  diagonal, 
this  yields. 


1 


A. -1 


1-A„  1  +  A(A.-1). 


^sphere  )  ' 


R,  =b 


(21) 


7?2  =  « (22-23) 


where  b  is  the  half-length  of  the  spheroid,  a  is  the  radius,  and  the  demagnetization  faetors  A„  are 
given  in  Smith  and  Morrison  (2006,  eqns.  A-2  and  A-3).  The  magneto-statie  response  is  simply 


= 


A.-l 


(24) 


l  +  4(//,-l) 

Combining  both  the  eleetro-dynamie  and  the  magneto-statie  responses,  we  obtain 


m: 


’  (0  =  2M„  (0  -  (t  =  5 . 5  5  •  1  O^O  -  <7 


(25) 


where  the  eleetro-dynamie  response  at  5.55  ms  is  subtraeted  to  aeeount  for  any  eddy  eurrents  that 
have  not  deeayed  when  the  transmitted  waveform  ehanges  slope.  Finally,  the  time  domain 
AFFTEM  response  is  ealeulated  from  Equations  3,  4,  and  25  where  co  is  replaeed  by  t. 
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The  inversion  employs  an  iterative  Gauss-Newton  minimization  eombined  with  step-size 
optimization  as  follows, 


aJ\x.)-  C  '  ■  (f  (x. )  -  /o )  =  x.+i  -  X. , 

(26) 

•^(a)  =  V/(x.), 

(27) 

(Tj  ...  0 

c= 

(28) 

- 1 

s: 

b 

o 

_ 1 

min/(4i)-/o".  (29) 

ae(0-l] 

where  is  the  pseudo-inverse  of  and  /(x)is  the  forward  model  veetor  eontaining 

predicted  data  samples,  /q  is  the  measured  data,  C  is  a  diagonal  matrix  containing  the  data 
standard  deviations  (it  is  assumed  that  the  data  are  identically  distributed  and  statistically 
independent),  and  x  contains  the  parameters.  At  each  iteration,  the  step  length  is  scaled  by  a 
which  is  varied  in  discrete  steps  over  the  interval  (0,1]  to  find  the  step  that  results  in  the  smallest 
least  squares  data  misfit.  Since  the  number  of  data  (number  of  coil  combinations  times  number 
of  instrument  locations  times  number  of  waveform  points)  far  exceeds  the  number  of  parameters 
to  be  estimated,  the  pseudo-inverse  of  the  Jacobian  matrix  is  calculated  by  singular  value 
decomposition.  Instabilities  due  to  poor  conditioning  are  avoided  by  scaling  the  parameters  so 
their  ranges  all  lie  within  1-2  orders  of  magnitude,  and  by  only  using  singular  values  that  have 
magnitudes  of  at  least  10'^  times  the  dominant  singular  value.  At  each  iteration,  the  target’s 
orientation  angles  are  constrained  to  be  between  -n  and  n,  and  the  spheroid  diameters  are 
constrained  to  be  positive.  Iteration  stops  when  a  local  minimum  is  found,  or  the  maximum 
allowable  number  of  iterations  is  reached.  This  basic  search  algorithm  is  used  with  a  state 
machine  to  estimate  the  target  parameters  as  described  below. 

The  first  step  in  the  parameter  estimation  problem  is  to  choose  an  initial  model.  The  x  and 
y  target  locations  are  selected  from  the  centroid  of  the  anomaly,  and  the  z  location  is  estimated 
by  finding  the  least  squares  minimum  of 


c  n 


r(p„,z)-r(o,z)1 
/;(p„,^)-/;(o,z)  j  ’ 


(30) 


where  /'’  is  the  modeled  data  and  is  the  measured  data  for  coil  combination  c  at  instrument 
location  n.  All  data  are  from  the  same  time-gate  (usually  t  =  5  ms),  and  pn  is  the  horizontal 
distance  from  the  centroid  at  instrument  location  n.  This  functional  calculates  the  slopes  of  the 
logarithm  of  the  fields  versus  distance  from  the  target  (i.e.,  slope  is  -3  for  the  fields  of  a  static 
dipole  decaying  at  f^),  and  compares  the  slopes  of  measured  data  and  data  modeled  for  a 


13 


spheroid.  To  determine  a  representative  permeability  value,  the  late-time  magnitudes  are 
examined  to  determine  if  the  target  is  ferrous  or  non-ferrous  (//^  =  !)•  For  ferrous  targets,  the 
demagnetization  effect  must  be  considered  (see  fig.  5).  For  targets  with  relative  permeabilities 
ranging  from  50  to  1000,  a  nominal  value  of  100  is  sufficient  to  model  objects  with  aspect  ratios 
less  than  about  four,  therefore  the  selected  representative  value  of  is  held  fixed  during  the 
inversion.  An  initial  conductivity  is  chosen  based  values  of  typical  metals  used  in  UXO 
construction.  The  conductivity  of  aluminum  alloys  typically  range  from  about  1.5-  lO’  to  3.5- 10^ 
S/m,  and  steel  alloys  typically  range  from  0.2-10’  to  0.9-10’  S/m,  which  makes  1.0- lO’  S/m  a 
reasonable  starting  value.  Initial  orientation  angles  (pitch  and  yaw)  are  zero. 

With  a  representative  permeability  value,  it  is  possible  to  determine  the  principle  spheroid 
diameters  using  the  response  at  time  t  =  0^  (the  instant  the  transmitter  turns  off).  This  however, 
requires  a  system  with  instantaneous  turn  off  time  and  a  receiver  with  infinite  bandwidth.  With 
some  high-bandwidth  systems,  it  may  be  possible  to  extrapolate  back  from  earliest  available  time 
sample  at  a  slope  of  t  to  estimate  the  dimensions  of  the  target. 

In  conducting  trials  with  the  minimization  algorithm,  it  was  observed  that  there  is  a  basin  of 
attraction  associated  with  both  prolate  models  and  oblate  models.  Oftentimes,  the  evolving  data 
misfit  function  would  enter  an  incorrect  basin  of  attraction  only  to  find  a  local  minimum.  The 
solution  to  this  problem  is  to  minimize  the  misfit  function  using  a  prolate  model,  then  using  an 
oblate  model,  and  then  choose  the  solution  with  the  best  fit.  While  minimizing  these  functions, 
only  data  from  a  single  time-gate  (typically  at  t  =  5  ms)  are  used,  and  both  initial  models  use  a 
larger  diameter  of  0.26  meter,  and  a  smaller  diameter  of  0.1  meter.  Because  most  of  the  energy 
in  the  ALLTEM  waveforms  from  ferrous  targets  at  times  greater  than  1  ms  is  magneto-static,  the 
fi  and  a  are  held  fixed  to  reduce  the  degrees  of  freedom  while  searching  for  optimum  prolate  and 
oblate  models.  The  final  step  is  to  polish  the  conductivity  value  using  the  best  solution  found 
thus  far,  using  data  from  all  time-gates  selected  for  analysis  and  holding  all  other  parameters 
fixed. 


The  mean  squared  error  in  the  best-fit  modeled  data  is  assumed  to  be  due  to  variations  from 
a  non-ideal  systematic  response.  These  variations  include  components  of  the  instrument 
response  not  accounted  for  by  the  model  (drift,  non-linear  response,  etc.),  components  of  target 
response  not  accounted  for  by  the  model,  ambient  EM  noise,  geologic  noise,  errors  in  instrument 
location,  and  attitude  variations  of  the  instrument.  To  estimate  the  uncertainty  in  the  estimated 
parameters,  each  parameter  is  perturbed  from  its  best-fit  value  until  the  mean  squared  error  of  the 
modeled  data  increases  by  the  variance  estimate  of  the  data. 

The  number  of  data  points  is  typically  chosen  to  be  less  than  -'1000  so  that  the  inversion 
can  be  accomplished  in  a  reasonable  time  frame  (about  a  minute).  When  selecting  a  set  of  coil 
combinations  to  use  in  the  analysis,  the  set  that  carries  the  most  (orthogonal)  information  is 
desirable.  To  select  a  subset  of  coil  combinations  from  the  recorded  set  of  19  coil  combinations, 
selections  are  made  in  order  of  decreasing  data  variance  until  a  single  selection  for  each  of  the 
nine  possible  polarization  combinations  (i.e.,  (Txx,  Rxx),  (Txx,  Rxy),  etc.)  has  been  made.  If 
more  coil  combinations  are  needed  to  fill  the  subset,  then  additional  selections  are  made  in  order 
of  decreasing  data  variance. 
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d.  What  tunable  parameters  (if  any)  are  used  in  the  characterization  process? 
(e.g.,  thresholds  on  background  noise,  etc.)?  Tunable  parameters  inelude  all  the  parameters 
derived  by  the  inversion  proeess. 

Classification:  This  seetion  should  include  the  details  describing  the  algorithm  and 
associated  data  and  parameters  used  for  discrimination  by  answering  the  following  questions: 

a.  What  algorithm  is  used  for  discrimination  (e.g.,  multi-layer  perception,  support  vector 
machine,  etc.)?  The  primary  algorithm  is  an  analysis  of  the  inverted  permeabilities  and  the 
lengths  and  widths  of  the  spheroid  models  produced  by  the  inversion  and  then  a  comparison  to 
coefficients  for  known  items  including  those  from  the  Calibration  grid  and  test  stand  data. 

b.  Why  is  this  algorithm  used  and  not  others?  This  discrimination  analysis  process  has 
been  used  successfully  for  the  ALLTEM  for  UXO  items. 

c.  Which  parameters  are  considered  as  possible  inputs  to  the  algorithm?  Permeabilities 
(composition  of  target  material)  and  dimensions  of  the  prolate  spheroid  resulting  from  the 
inversion  analysis. 

d.  What  are  the  outputs  of  the  algorithm  (probabilities,  confidence  levels)?  Multiple 
probabilities  of  classification  with  associated  confidence  levels  are  derived  for  a  given  target 
item.  These  probabilities  represent  the  likelihood  of  an  item  being  clutter  or  ordnance  and  the 
most  likely  types  of  ordnance. 

e.  How  is  the  threshold  set  to  decide  where  the  munitions /non-munitions  line  lies  in  the 
discrimination  process?  The  threshold  used  to  determine  UXO  vs  clutter  is  based  on  the  ratio  of 
the  spheroid  dimensions.  For  a  rod-like  item,  the  two  smaller  widths  should  be  similar  and  much 
smaller  than  the  third,  much  larger,  length.  Clutter  typically  does  not  follow  this  pattern  although 
some  can. 

Training:  This  section  should  include  the  details  of  how  training  data  is  used  to  make  a 
decision  on  the  likelihood  of  the  anomaly  correspondence  to  munitions.  Please  answer  the 
following  questions: 

a.  Which  tunable  parameters  have  final  values  that  are  optimized  over  a  training  set  of 
data  and  which  have  values  that  are  set  according  to  geophysical  knowledge  (i.e.,  intuition, 
experience,  common  sense)?  Training  data  is  used  to  tune  estimates  of  location,  depth,  spheroid 
dimensions,  time  decay  constants,  and  composition  analysis.  Geophysical  knowledge  comes  in 
when  deciding  that  a  rod-like,  sphere-like,  or  disk-like  object  is  a  UXO  versus  a  piece  of  clutter. 

(1)  For  those  tunable  parameters  with  final  values  set  according  to  geophysical 
knowledge: 

(a)  What  is  the  reasoning  behind  choosing  these  particular  values?  These  shapes  (rod, 
sphere,  disk)  seem  to  be  the  typical  type  of  ordnance  used  on  training  ranges. 
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(b)  Why  were  the  final  values  not  optimized  over  a  training  set  of  data?  While  the  final 
values  are  optimized  over  a  training  set  of  data,  they  are,  to  a  large  degree,  based  on  a  priori  site 
specific  data  at  a  given  site. 

(2)  For  those  tunable  parameters  with  final  values  optimized  over  the  training  set  data: 

(a)  What  training  data  is  used  (e.g.,  all  data,  a  randomly  chosen  portion  of  data,  etc.)? 
All  available  data  is  utilized  to  train  the  inversion  and  classification  algorithms. 

(b)  What  error  metric  is  minimized  during  training  (e.g.,  mean  squared  error,  etc.)? 
Inferred  composition  analysis  and  definition  of  an  ordnance  by  its  spheroid  equivalent 
dimensions  and  time  decay.. 

(c)  What  learning  rule  is  used  during  training  (e.g.,  gradient  descent,  etc.)?  Determine 
best  parameters  to  identify  and  characterize  various  ordnance  versus  clutter. 

(d)  What  criterion  is  used  to  stop  training  (e.g.,  number  of  iterations  exceeds  threshold, 
good  generalization  over  validation  set  of  data,  etc.)?  Criterion  is  limit  of  the  number  of  training 
items.  The  more  known  items  available  the  better  the  statistical  analysis  of  the  possible 
variations  for  a  given  type  of  ordnance. 

(e)  Are  all  tunable  parameters  optimized  at  once  or  in  sequence 
(“in  sequence  ”  =  parameters  1  is  held  constant  at  some  common  sense  values  while  parameter  2 
is  optimized,  and  then  parameter  2  is  held  constant  at  its  optimized  value  while  parameter  1  is 
optimized)?  Tunable  parameters  are  optimized  in  sequence. 

b.  What  are  the  final  values  of  all  tunable  parameters  for  the  characterization  process? 
The  final  values  for  the  characterization  are  the  correctly  classified  targets. 

2,1,4  Data  Submission  Format 


Data  were  submitted  for  scoring  in  accordance  with  data  submission  protocols  outlined  on 
the  USAEC  Web  site  www.uxotestsites.org.  These  submitted  data  are  not  included  in  this  report 
in  order  to  protect  GT  information. 

2,1,5  Demonstrator  Quality  Assurance  (QA)  and  Quality  Control  IOC)  (provided  by 

demonstrator) 

Overview  of  Quality  Control  (QC):  The  ALLTEM  system  has  a  real-time  data  display  that 
instantly  shows  the  operator  if  the  transmitting/receiving  functions  of  the  system  fail.  In 
addition,  we  plan  to  find  a  location  with  no  known  targets  and  repetitively  reoccupy  that  location 
and  record  data,  including  GPS  data,  to  assess  and  document  any  drifts  that  may  occur  in  the 
instrumentation.  Standard  operating  procedure  with  all  these  systems  is  to  occupy  a  designated 
clean  location  at  least  twice  each  day:  prior  to  and  at  the  completion  of  regular  data  acquisition. 
This  usually  takes  place  in  the  morning  and  afternoon,  but  in  case  of  an  extended  pause  in  the 
middle  of  the  day,  an  additional  reference  data  set  may  be  acquired.  This  will  also  test  the 
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accuracy  and  repeatability  of  the  navigation  data.  As  with  all  analog  and  time-base  systems,  drift 
will  occur  mainly  due  to  component  tolerances  and  temperature  dependencies.  This  inherent 
system  drift  limits  the  absolute  accuracy  of  the  measurements  that  can  be  attained.  The  reference 
data  are  used  primarily  as  a  metrie  for  overall  aeeuraey.  Abnormal  drift,  as  would  be  eaused  by 
battery  depletion  or  component  degradation,  would  trigger  a  system  eheck  and  data  review.  The 
hardware  problem  would  be  corrected  and  field  data  acquisition  would  resume.  Any  previous 
data  deemed  degraded  would  be  reacquired.  We  also  plan  to  preprocess  data  overnight  or 
eoneurrent  with  data  aequisition  to  visually  ensure  that  there  are  no  serious  “glitehes”  or  “tears” 
in  the  data.  Any  corrupted  lines  will  be  repeated.  The  GPS  will  be  referenced  to  a  loeal  geodetic 
marker. 

Overview  of  Quality  Assurance  (QA):  As  mentioned  above,  the  planned  along-line  data 
density  will  be  around  15  to  20  cm  with  a  line  spacing  of  50  cm.  This  will  ensure  that  the 
1 -meter  square  antennas  will  sample  over  every  point  on  the  ground.  The  basie  position  aeouraey 
of  our  real-time  kinematic  differential  GPS  system  is  better  than  2  em  when  operating  in  “fixed” 
mode.  The  Lab  VIEW  program  reads  the  GPS  data  and  mode.  If  the  mode  is  not  fixed,  the 
Lab  VIEW  program  flashes  a  visual  warning  on  the  monitor  to  alert  the  operator  that  the  GPS  is 
not  in  fixed  mode.  However,  after  our  reeent  experienee  with  radio  eommunications  problems  at 
YPG,  we  now  reeord  to  memory  eards  the  raw  data  at  both  the  GPS  rover  and  base  stations  and 
then  post-process  the  data  to  correct  locations.  Other  sources  of  error  in  positioning,  such  as  GPS 
data  latency,  GPS  antenna-to-sensor  offset,  and  tilting  of  the  GPS  antenna  mast  with  topography 
degrade  absolute  position  aeeuraey.  We  have  added  an  Attitude  Heading  and  Reference  System 
(AHRS)  to  measure  the  cart  orientation  relative  to  the  ground.  We  have  also  developed  a 
navigation  program  in  LabView  that  runs  concurrent  with  the  aequisition  program  to  maintain 
position  over  large  distances. 

Data  proeessing  will  begin  in  the  field.  At  the  end  of  each  survey  line,  the  data  is  automatically 
copied  to  an  external  hard  drive  which  will  be  swapped  out  with  another  drive  periodically 
during  the  survey.  The  data  is  then  quickly  batch  processed  in  Geosoft  Oasis  Montaj  and  within 
minutes  the  quality  of  the  survey  data  density  and  areal  eoverage  ean  be  evaluated. 

2,1,6  Additional  Records 

The  following  reoord(s)  by  this  vendor  can  be  aceessed  via  the  Internet  as  Mierosoft  Word 
doeuments  at  www.uxotestsites.org. 
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2.2  APG  SITE  INFORMATION 


2,2,1  Location 


The  APG  Standardized  Test  Site  is  located  within  a  secured  range  area  of  the  Aberdeen 
Area.  The  Aberdeen  Area  of  APG  is  located  approximately  30  miles  northeast  of  Baltimore  at 
the  northern  end  of  the  Chesapeake  Bay.  The  Standardized  Test  Site  encompasses  17  acres  of 
upland  and  lowland  flats,  woods,  and  wetlands. 

2.2.2  Soil  Type 

According  to  the  soils  survey  conducted  for  the  entire  area  of  APG  in  1998,  the  test  site 
consists  primarily  of  Elkton  Series  type  soil  (ref  2).  The  Elkton  Series  consist  of  very  deep, 
slowly  permeable,  poorly  drained  soils.  These  soils  formed  in  silty  aeolin  sediments  and  the 
underlying  loamy  alluvial  and  marine  sediments.  They  are  on  upland  and  lowland  flats  and  in 
depressions  of  the  Mid-Atlantic  Coastal  Plain.  Slopes  range  from  0  to  2  percent. 

ERDC  conducted  a  site-specific  analysis  in  May  2002  (ref  3).  The  results  basically 
matched  the  soil  survey  mentioned  above.  Seventy  percent  of  the  samples  taken  were  classified 
as  silty  loam.  The  majority  (77  percent)  of  the  soil  samples  had  a  measured  water  content 
between  15  and  30  percent  with  the  water  content  decreasing  slightly  with  depth. 

Eor  more  details  concerning  the  soil  properties  at  the  APG  test  site,  go  to 
www.uxotestsites.org  on  the  Web  to  view  the  entire  soils  description  report. 

2.2.3  Test  Areas 


A  description  of  the  test  site  areas  at  APG  is  presented  in  Table  1.  A  test  site  layout  is 
shown  in  Eigure  2. 
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TABLE  1,  TEST  SITE  AREAS 


Area 

Description 

Calibration  lanes 

Contains  14  standard  munitions  items  buried  in  six  positions,  with  representation 
of  clutter,  at  various  angles  and  depths  to  allow  demonstrators  to  calibrate  their 
equipment. 

Blind  grid 

Contains  400  grid  cells  in  a  0.5-acre  site.  The  center  of  each  grid  cell  contains 
either  munitions,  clutter,  or  nothing. 

Open  field 

A  10-acre  site  composed  of  generally  open  and  flat  terrain  with  minimal  clutter 
and  minor  navigational  obstacles.  Vegetation  height  varies  from  15  to  25  cm. 
This  area  is  subdivided  into  four  subareas  (legacy,  direct  fire,  indirect  fire,  and 
challenge). 

•  Open  field  (legacy) 

The  legacy  subarea  contains  the  same  wide  variety  of  randomly-placed  munitions 
that  were  present  in  the  open  field  prior  to  the  January  2008  general 
reconfiguration  of  the  site. 

•  Open  field  (direct  fire) 

The  direct  fire  subarea  contains  only  three  munition  types  that  could  be  typically 
found  at  an  impact  area  of  a  direct  fire  weapons  range.  Munitions  and  clutter  are 
placed  in  a  pattern  typical  for  these  munitions. 

•  Open  field  (indirect  fire) 

The  indirect  fire  subarea  contains  only  three  munition  types  that  could  be  typically 
found  at  an  impact  area  of  an  indirect  fire  weapons  range.  Munitions  and  clutter 
are  placed  in  a  pattern  typical  for  these  munitions. 

•  Open  field  (challenge) 

The  challenge  subarea  is  easily  reconfigurable  to  meet  the  specific  needs  and 
requirements  of  the  demonstrator  or  the  program  sponsor.  Any  results  from  this 
area  are  not  reported  in  the  standardized  scoring  record. 

Woods 

1.34-acre  area  consisting  of  cleared  woods  (tree  removal  with  only  stumps 
remaining),  partially  cleared  woods  (including  all  underbrush  and  fallen  trees), 
and  virgin  woods  (i.e.,  woods  in  natural  state  with  all  trees,  underbrush,  and 
fallen  trees  left  in  place). 

Moguls 

1.30-acre  area  consisting  of  two  areas  (the  rectangular  or  driving  portion  of  the 
course  and  the  triangular  section  with  more  difficult,  nondrivable  terrain).  A 
series  of  craters  (as  deep  as  0.91  m)  and  mounds  (as  high  as  0.91  m)  encompass 
this  section. 
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Figure  2.  Test  site  layout. 


20 


2,2,4  Standard  and  Nonstandard  Inert  Munitions  Targets 

The  standard  and  nonstandard  munitions  items  emplaeed  in  the  test  areas  are  presented  in 
Table  2.  Standardized  targets  are  members  of  a  set  of  speeific  munitions  items  that  have 
identieal  properties  to  all  other  items  in  the  set  (ealiber,  eonfiguration,  size,  weight,  aspeet  ratio, 
material,  filler,  magnetie  remanenee,  and  nomenelature).  Nonstandard  targets  are  inert 
munitions  items  having  properties  that  differ  from  those  in  the  set  of  standardized  items. 
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TABLE  2,  INERT  MUNITIONS  TARGETS 


Munition 

Calibration 

Open  Field 

Open  Field 

Open  Field 

Item 

Type 

Lanes 

Blind  Grid 

Direct  Fire 

Indirect  Fire 

Legacy 

Moguls 

Woods 

20-mm  Projectile  M55 

S 

X 

X 

X 

X 

25-mm  Projectile  M794 

S 

X 

X 

X 

3 7 -mm  Projectile  M47 

s 

X 

X 

X 

40-mm  Projectile  MKII  Bodies 

s 

X 

X 

X 

X 

BDU-28  Submunition 

s 

X 

X 

X 

X 

BLU-26  Submunition 

s 

X 

X 

X 

X 

M42  Submunition 

s 

X 

X 

X 

X 

57-mm  Projectile  APC  M86 

s 

X 

X 

X 

X 

60-mm  Mortar  M49A3 

s 

X 

X 

X 

2.75-in.  Rocket  M230 

s 

X 

X 

X 

X 

81 -mm  Mortar  M3  74 

s 

X 

X 

X 

X 

X 

X 

105-mm  HEAT  Rounds  M456 

s 

X 

X 

X 

105 -mm  HEAT  Round  M490 

s 

X 

X 

X 

105-mm  Projectile  M60 

s 

X 

X 

X 

X 

X 

X 

155-mm  Projectile  M483A1 

s 

X 

X 

X 

X 

20-mm  Projectile  M55 

NS 

X 

X 

X 

20-mm  Projectile  M97 

NS 

X 

X 

X 

40-mm  Projectile  M813 

NS 

X 

X 

X 

60-mm  Mortar  (JPG) 

NS 

X 

X 

X 

60-mm  Mortar  M49 

NS 

X 

X 

X 

2.75-in.  Rocket  M230 

NS 

X 

X 

X 

2.75-in.  Rocket  XM229 

NS 

X 

X 

X 

81 -mm  Mortar  (JPG) 

NS 

X 

X 

X 

8 1  -mm  Mortar  M3  74 

NS 

X 

X 

X 

105-mm  Projectile  M60 

NS 

X 

X 

X 

155-mm  Projectile  M483A 

NS 

X 

X 

X 

HEAT  =  high-explosive  antitank. 

JPG  =  Jefferson  Proving  Ground. 

NS  =  Nonstandard  munition. 

S  =  Standard  munition. 


2,3  ATC  SURVEY  COMMENTS 


Survey  of  the  Open  Field  (Indireet  Fire  Area)  was  only  partially  eompleted. 
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SECTIONS.  FIELD  DATA 


3.1  DATE  OF  FIELD  ACTIVITIES  (2  to5,  8  tolO  March  2010) 

3.2  AREAS  TESTED/NUMBER  OF  HOURS 

Areas  tested  and  total  numbers  of  hours  operated  at  eaeh  site  are  presented  in  Table  3. 


TABLES.  AREAS  TESTED  AND 
NUMBER  OF  HOURS 


Area 

Number  of  Hours 

Calibration  lanes 

4.08 

Blind  grid 

5.92 

Open  field 

36.33 

Woods 

- 

Mogul 

- 

Mine  grid 

- 

Note;  Table  3  represents  the  total  time  spent  in  eaeh  area. 


3.3  TEST  CONDITIONS 
3.3.1  Weather  Conditions 

An  APG  weather  station  located  approximately  1  mile  west  of  the  test  site  was  used  to 
record  average  temperature  and  preeipitation  on  a  half-hour  basis  for  each  day  of  operation.  The 
temperatures  presented  in  Table  4  represent  the  average  temperature  during  field  operations  from 
0700  to  1700  hours,  while  preeipitation  data  represent  a  daily  total  amount  of  rainfall.  Hourly 
weather  logs  used  to  generate  this  summary  are  provided  in  Appendix  B. 


TABLE  4.  TEMPERATURE/PRECIPITATION  DATA  SUMMARY 


Date,  2010 

Average  Temperature,  “F 

Total  Daily  Precipitation,  in. 

2  March 

40.4 

0.08 

3  March 

40.4 

0.02 

4  March 

41.8 

0.00 

5  March 

43.1 

0.00 

8  March 

51.8 

0.00 

9  March 

55.0 

0.00 

1 0  March 

50.9 

0.00 
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3,3,2  Field  Conditions 


USGS  surveyed  the  calibration  grid,  blind  grid,  direct  and  indirect  fire  areas.  The  field  was 
wet  in  all  areas  of  the  field  due  to  rain  and  snow  prior  to  testing.  USGS  was  unable  to  complete  a 
small  portion  of  the  indirect  fire  area  due  to  standing  water. 

3,3,3  Soil  Moisture 


Three  soil  probes  were  placed  at  various  locations  within  the  site  to  capture  soil  moisture 
data:  blind  grid,  calibration,  open  field,  and  wooded  areas.  Measurements  were  collected  in 
percent  moisture  and  were  taken  twice  daily  (morning  and  afternoon)  from  five  different  soil 
depths  (1  to  6  in.,  6  to  12  in.,  12  to  24  in.,  24  to  36  in.,  and  36  to  48  in.)  from  each  probe.  Soil 
moisture  logs  are  provided  in  Appendix  C. 

3,4  FIELD  ACTIVITIES 

3.4.1  Setup/Mobilization 

These  activities  included  initial  mobilization  and  daily  equipment  preparation  and 
breakdown.  A  three-person  crew  took  5  hours  and  40  minutes  to  perform  the  initial  setup  and 
mobilization.  A  total  of  5  hours  and  20  minutes  of  equipment  preparation  was  accrued,  and  end 
of  day  equipment  breakdown  totaled  2  hours. 

3.4.2  Calibration 


USGS  spent  a  total  of  4  hours  5  minutes  in  the  calibration  lanes,  of  which  2  hours  and 
30  minutes  were  spent  collecting  data.  Numerous  calibration  exercises  occurred  while  surveying 
the  Blind  Grid  and  Open  Field  lasting  4  hours  and  5  minutes. 

3,4,3  Downtime  Occasions 


Occasions  of  downtime  are  grouped  into  five  categories:  equipment/data  checks  or 
equipment  maintenance,  equipment  failure  and  repair,  weather,  demonstration  site  issues,  or 
breaks/lunch.  All  downtime  is  included  for  the  purposes  of  calculating  labor  requirements 
(section  5)  except  for  downtime  due  to  demonstration  site  issues.  Demonstration  site  issues, 
while  noted  in  the  daily  log,  are  considered  nonchargeable  downtime  for  the  purposes  of 
calculating  labor  costs  and  are  not  discussed.  Breaks  and  lunches  are  discussed  in  this  section 
and  billed  to  the  total  site  survey  area. 

3.4.3. 1  Equipment/data  checks,  maintenance.  Equipment  data  checks  and  maintenance 
activities  accounted  for  1  hour  and  35  minutes  of  site  usage  time.  These  activities  included 
changing  out  batteries  and  performing  routine  data  checks  to  ensure  the  data  were  being  properly 
recorded/collected.  USGS  spent  1  hour  and  20  minutes  for  breaks  and  lunches. 

3.4.3.2  Equipment  failure  or  repair.  Two  equipment  failures  occurred  during  this  survey 
totaling  1  hour  and  40  minutes.  A  faulty  GPS  cable  had  to  be  repaired.  Also,  USGS  had  two 
incidents  of  flat  tires  while  surveying.  All  repairs  were  made  on  site. 

3.4.3.3  Weather,  No  weather  delays  occurred  during  the  survey. 
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3,4,4  Data  Collection 


TABLES,  TOTAL  TIME 
NRL,  SPENT  PER  AREA 


Area 

Time,  hr/min 

Blind  grid 

3  hours/50  minutes 

Open  field 

— 

Legacy 

— 

Direct  fire 

8  hours/ 15  minutes 

Indirect  fire 

13  hours/5  minutes 

Challenge 

— 

Wooded 

— 

Mine  Grid 

— 

Moguls 

— 

Note;  Table  5  represents  the  total  time  spent  in  eaeh  area  colleeting  data. 


3,4,5  Demobilization 


The  uses  survey  erew  eondueted  a  demonstration  of  the  ealibration,  blind  grid,  and  open 
field,  indirect/direet  fire  areas.  Demobilization  oecurred  on  10  March  2010.  On  that  day,  it  took 
the  crew  3  hours  and  10  minutes  to  break  down  and  pack  up  their  equipment. 

3.5  PROCESSING  TIME 

USGS  submitted  the  raw  data  from  the  demonstration  activities  on  the  last  day  of  the 
demonstration,  as  required.  The  initial  scoring  submittal  data  was  provided  September  2010. 
However,  the  final  submittal  for  scoring  was  provided  March  2011. 

3.6  DEMONSTRATOR’S  FIELD  PERSONNEL 

Ted  Asch 
Craig  Moulton 
Moustapha  Sylla 

3.7  DEMONSTRATOR’S  FIELD  SURVEYING  METHOD 

USGS  collected  the  data  in  a  linear  fashion,  using  a  line  spacing  the  width  of  the  array. 

3.8  SUMMARY  OF  DAILY  LOGS 

Daily  logs  capture  all  field  activities  during  this  demonstration  and  are  provided  in 
Appendix  D. 
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SECTION  4.  TECHNICAL  PERFORMANCE  RESULTS 


4.1  ROC  CURVES  USING  ALL  MUNITIONS  CATEGORIES 

The  probability  of  detection  for  the  response  stage  and  the  discrimination  stage 

(p^disc)  respective  probability  of  clutter  detection  or  probability  of  false  positive 

within  each  area  are  shown  in  Figures  3  through  8.  The  probabilities  plotted  against 
their  respective  background  alarm  rate  within  each  area  are  shown  in  Figures  9  through  14. 
Both  figures  use  horizontal  lines  to  illustrate  the  performance  of  the  demonstrator  at  two 
demonstrator-specified  points;  at  the  system  noise  level  for  the  response  stage,  representing  the 
point  below  which  targets  are  not  considered  detectable,  and  at  the  demonstrator’s  recommended 
threshold  level  for  the  discrimination  stage,  defining  the  subset  of  targets  the  demonstrator  would 
recommend  digging  based  on  discrimination. 


-  Response 

- Discrimination 

---  Noise  Level 
^—Threshold 


Figure  3.  ALLTEM/towed  blind  grid  probability  of  detection  for  response  and  discrimination 
stages  versus  their  respective  probability  of  false  positive. 
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-  Response 

- Discrimination 

---  Noise  Level 
- Threshold 


Figure  4.  ALLTEM/towed  open  field  (direct  fire)  probability  of  detection  for  response  and 
discrimination  stages  versus  their  respective  probability  of  false  positive. 


-  Response 

- Discrimination 

- Noise  Level 

- Threshold 


Figure  5.  ALLTEM/towed  open  field  (indirect  fire)  probability  of  detection  for  response  and 
discrimination  stages  versus  their  respective  probability  of  false  positive. 
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Not  covered 


Figure  6.  ALLTEM/towed  open  field  (legacy)  probability  of  detection  for  response  and 
diserimination  stages  versus  their  respeetive  probability  of  false  positive. 


Not  covered 

Figure  7.  ALLTEM/towed  wooded  probability  of  deteetion  for  response  and  discrimination 
stages  versus  their  respective  probability  of  false  positive. 


Not  covered 

Figure  8.  ALLTEM/towed  mogul  probability  of  deteetion  for  response  and  diserimination 
stages  versus  their  respective  probability  of  false  positive. 


-  Response 

- Discrimination 

- Noise  Level 

- Threshold 


Figure  9.  ALLTEM/towed  blind  grid  probability  of  detection  for  response  and  discrimination 
stages  versus  their  respeetive  probability  of  baekground  alarm. 
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- Response 

- Discrimination 

- Noise  Level 

- Threshold 


Background  Alarm  Rate 


Figure  10.  ALLTEM/towed  open  field  (direct  fire)  probability  of  detection  for  response  and 
discrimination  stages  versus  their  respective  background  alarm  rate. 


- Response 

- Discrimination 

- Noise  Level 

- Threshold 


Figure  1 1 .  ALLTEM/towed  open  field  (indirect  fire)  probability  of  detection  for  response  and 
discrimination  stages  versus  their  respective  background  alarm  rate. 
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Not  covered 


Figure  12.  ALLTEM/towed  open  field  (legacy)  probability  of  detection  for  response  and 
diserimination  stages  versus  their  respeetive  baekground  alarm  rate. 

Not  covered 

Figure  13.  ALLTEM/towed  wooded  probability  of  deteetion  for  response  and  diserimination 
stages  versus  their  respective  background  alarm  rate. 

Not  eovered 

Figure  14.  ALLTEM/towed  mogul  probability  of  detection  for  response  and  discrimination 
stages  versus  their  respective  baekground  alarm  rate. 


4.2  PERFORMANCE  SUMMARIES 

Results  for  each  of  the  testing  areas  are  presented  in  Tables  6  (for  labor  requirements,  see 
section  5).  The  response  stage  results  are  derived  from  the  list  of  anomalies  above  the 
demonstrator-provided  noise  level.  The  results  for  the  diserimination  stage  are  derived  from  the 
demonstrator’s  reeommended  threshold  for  optimizing  munitions  related  cleanup  by  minimizing 
false  alarm  digs  and  maximizing  munitions  recovery.  The  lower  and  upper  90-pereent 
confidence  limits  on  Pd,  Pcd,  and  Pfj,  were  calculated  assuming  that  the  number  of  deteetions  and 
false  positives  are  binomially  distributed  random  variables. 
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TABLE  6a.  BLIND  GRID  TEST  AREA  RESULTS 


1  Response  Stage 

Discrimination  Stage 

Munitions’* 

Scores 

P/“:  by  type 

All  Types 

105-mm 

81/60-mm 

37/25-mm 

All  Types 

105-mm 

81/60-mm 

37/25-mm 

1.00 

1.00 

0.98 

1.00 

1.00 

0.93 

1.00 

1.00 

0.93 

1.00 

1.00 

0.93 

0.98 

0.96 

0.91 

1.00 

0.97 

0.88 

0.98 

0.93 

0.83 

1.00 

0.97 

0.88 

1  By  Depth^  \ 

0  to  4D 

1.00 

1.00 

1.00 

1.00 

0.96 

1.00 

0.92 

1.00 

4D  to  8D 

1.00 

1.00 

1.00 

1.00 

0.97 

1.00 

1.00 

0.95 

8D  to  12D 

1.00 

1.00 

1.00 

1.00 

0.89 

0.83 

1.00 

1.00 

Clutter 

Scores 

Pci 

PfP 

1  By  Mass  \ 

By  Depth'’ 

All  Mass 

0  to  0.25  kg 

>0.25  to 

1  kg 

>1  to  8  kg 

All  Mass 

0  to  0.25  kg 

>0.25  to 

1  kg 

>1  to  8  kg 

All  Depth 

0.92 

0.88 

0.84 

0.76 

1.00 

1.00 

0.27 

0.22 

0.17 

0.12 

0.27 

0.50 

0  to  0.15  m 

0.89 

0.77 

1.00 

1.00 

0.21 

0.11 

0.27 

0.67 

0.15  to  03  m 

0.88 

0.60 

1.00 

1.00 

0.25 

0.20 

0.29 

0.25 

0.3  to  0.6  m 

NA 

NA 

NA 

NA 

NA 

NA 

NA 

NA 

Background  Alarm  Rates 

Pb/“:  0.24  1  0.04 

^In  cells  with  offset  data  entries,  the  numbers  to  the  left  are  the  result  and  the  two  numbers  to  the 
right  are  an  upper  and  lower  90-pereent  confidenee  interval  for  an  assumed  binomial 
distribution. 

’’All  depths  are  measured  to  the  eenter  of  the  objeet. 
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TABLE  6b.  OPEN  FIELD  DIRECT  FIRE  TEST  AREA  RESULTS 


1  Response  Stage 

Discrimination  Stage 

Munitions'' 

Scores 

P/“.-  by  type 

All  Types 

105-mm 

37-mm 

25-mm 

All  Types 

105-mm 

37-mm 

25-mm 

0.97 

0.95 

0.93 

0.98 

0.95 

0.91 

0.97 

0.93 

0.89 

0.99 

0.97 

0.93 

0.95 

0.93 

0.90 

0.95 

0.92 

0.86 

0.96 

0.92 

0.87 

0.97 

0.94 

0.89 

1  By  Density  \ 

High 

0.94 

0.96 

0.89 

0.96 

0.91 

0.96 

0.85 

0.92 

Medium 

0.97 

1.00 

0.93 

0.97 

0.93 

0.93 

0.93 

0.94 

Low 

0.95 

0.90 

0.97 

0.97 

0.93 

0.86 

0.97 

0.94 

1  By  Depth^  \ 

0  to  4D 

0.96 

1.00 

0.94 

1.00 

0.94 

1.00 

0.92 

0.92 

4D  to  8D 

0.96 

0.96 

0.92 

0.98 

0.94 

0.93 

0.92 

0.96 

8D  to  12D 

0.91 

0.89 

0.00 

0.92 

0.84 

0.78 

0.00 

0.88 

Clutter 

Scores 

Pci 

PfP 

1  By  Mass  \ 

By  Depth'’ 

All  Mass 

0  to  0.25 1^ 

>0.25  to 

1  kg 

>1  to  8  kg 

All  Mass 

0  to  0.25 1^ 

>0.25  to 

1  kg 

>1  to  8  kg 

All  Depth 

0.76 

0.73 

0.69 

0.57 

0.86 

0.96 

0.36 

0.32 

0.28 

0.13 

0.46 

0.52 

0  to  0.15  m 

0.73 

0.59 

0.87 

0.95 

0.31 

0.14 

0.50 

0.45 

0.15  to  03  m 

0.69 

0.54 

0.75 

1.00 

0.24 

0.08 

0.13 

0.60 

0.3  to  0.6  m 

0.80 

0.00 

1.00 

1.00 

0.70 

0.00 

0.00 

1.00 

Background  Alarm  Rates 

BAR™*:  0.11  1  BAR"'™:  0.04 

Groups 

Found 

0.94 

0.88 

Identifled 

0.12 

0.06 

Coverage 

0.53 

0.47 

^In  cells  with  offset  data  entries,  the  numbers  to  the  left  are  the  result  and  the  two  numbers  to  the 
right  are  an  upper  and  lower  90-pereent  eonfidenee  interval  for  an  assumed  binomial 
distribution. 

’’All  depths  are  measured  to  the  eenter  of  the  objeet. 
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TABLE  6c.  OPEN  FIELD  INDIRECT  FIRE  TEST  AREA  RESULTS 


1  Response  Stage 

Discrimination  Stage 

Munitions'* 

Scores 

1  pr:  by  type  1 

by  type 

All  Types 

105-mm 

81-mm 

60-mm 

All  Types 

105-mm 

81-mm 

60-mm 

0.96 

0.95 

0.92 

0.99 

0.97 

0.92 

0.96 

0.93 

0.87 

0.97 

0.94 

0.90 

0.92 

0.89 

0.86 

0.97 

0.94 

0.90 

0.87 

0.82 

0.75 

0.95 

0.91 

0.86 

1  By  Density  \ 

High 

0.90 

0.92 

0.86 

0.92 

0.86 

0.92 

0.86 

0.80 

Medium 

0.93 

1.00 

0.89 

0.90 

0.87 

0.97 

0.75 

0.90 

Low 

0.99 

0.97 

1.00 

1.00 

0.93 

0.93 

0.84 

1.00 

1  By  Depth^  \ 

0  to  4D 

0.97 

1.00 

0.97 

0.95 

0.93 

0.96 

0.88 

0.92 

4D  to  8D 

0.93 

0.97 

0.90 

0.93 

0.91 

0.97 

0.85 

0.93 

8D  to  12D 

0.84 

0.50 

0.89 

0.92 

0.64 

0.50 

0.44 

0.83 

Clutter 

Scores 

PfP 

1  By  Mass  \ 

By  Depth'’ 

All  Mass 

0  to  0.25 1^ 

>0.25  to 

1  kg 

>1  to  8  kg 

All  Mass 

0  to  0.25  kg 

>0.25  to 

1  kg 

>1  to  8  kg 

All  Depth 

0.68 

0.64 

0.61 

0.50 

0.79 

0.84 

0.08 

0.06 

0.05 

0.04 

0.05 

0.21 

0  to  0.15  m 

0.62 

0.50 

0.78 

0.83 

0.07 

0.05 

0.06 

0.35 

0.15  to  03  m 

0.75 

0.63 

0.79 

0.82 

0.02 

0.00 

0.00 

0.06 

0.3  to  0.6  m 

0.82 

0.00 

0.83 

1.00 

0.00 

0.00 

0.00 

0.00 

Background  Alarm  Rates 

BAR™*:  0.27  |  BAR******:  0.02 

Groups 

Found 

0.95 

0.80 

Identified 

0.00 

0.00 

Coverage 

0.47 

0.40 

^In  cells  with  offset  data  entries,  the  numbers  to  the  left  are  the  result  and  the  two  numbers  to  the 
right  are  an  upper  and  lower  90-pereent  eonfidenee  interval  for  an  assumed  binomial 
distribution. 

’’All  depths  are  measured  to  the  eenter  of  the  objeet. 


Note;  This  Results  Table  was  seored  using  a  subset  of  the  full  Open  Field  (Indireet  Fire) 
Ground  Truth,  therefore  these  results  eannot  be  direetly  eompared  to  any  other  APG 
Open  Field  (Indireet  Fire)  Results  reported. 
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TABLE  6d.  OPEN  FIELD  LEGACY  TEST  AREA  RESULTS  (not  covered) 


1  Response  Stage 

Discrimination  Stage 

Munitions’* 

Scores 

by  type 

All  Types 

Small 

Medium 

Large 

All  Types 

Small 

Medium 

Large 

- 

- 

-- 

- 

- 

- 

- 

- 

1  By  Depth^  \ 

0  to  4D 

- 

- 

- 

- 

- 

- 

- 

- 

4D  to  8D 

- 

- 

- 

- 

- 

- 

- 

- 

8D  to  12D 

- 

- 

- 

- 

- 

- 

- 

- 

>  12D 

- 

- 

- 

- 

- 

- 

- 

- 

Clutter 

Scores 

Pci 

PfP 

1  By  Mass  \ 

By  Depth^ 

All 

Mass 

Oto 

0.25  kg 

>0.25  to 

1  kg 

>1  to 

10  kg 

>  10  kg 

All 

Mass 

Oto 

0.25  kg 

>0.25  to 

1  kg 

>1  to 

8  kg 

<10kg 

All  Depth 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

0  to  0.15  m 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

0.15to0.3m 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

0.3  to  0.6  m 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

>  0.6  m 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

Background  Alarm  Rates 

BAR™*:  1  BAR""**^: 

Groups 

Found 

- 

- 

Identifled 

- 

- 

Coverage 

- 

- 

^In  cells  with  offset  data  entries,  the  numbers  to  the  left  are  the  result  and  the  two  numbers  to  the 
right  are  an  upper  and  lower  90-pereent  eonfidence  interval  for  an  assumed  binomial 
distribution. 

’’All  depths  are  measured  to  the  center  of  the  object. 
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TABLE  6e.  WOODED  TEST  AREA  RESULTS  (not  covered) 


1  Response  Stage 

Discrimination  Stage 

Munitions’* 

Scores 

P/'‘;  by  type 

All  Types 

Small 

Medium 

Large 

All  Types 

Small 

Medium 

Large 

- 

- 

- 

- 

- 

- 

- 

- 

1  By  Depth^  \ 

0  to  4D 

- 

- 

- 

- 

- 

- 

- 

- 

4D  to  8D 

- 

- 

- 

- 

- 

- 

- 

- 

8D  to  12D 

- 

- 

- 

- 

- 

- 

- 

- 

>  12D 

- 

- 

- 

- 

- 

- 

- 

- 

Clutter 

Scores 

Pci 

PfP 

1  By  Mass  \ 

All 

Mass 

Oto 

0.25  kg 

>0.25  to 

1  kg 

>1  to 

10  kg 

>10  kg 

All 

Mass 

Oto 

0.25  kg 

>0.25  to 

1  kg 

>1  to 

8  kg 

<  10kg 

All  Depth 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

0  to  0.15  m 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

0.15to0.3m 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

0.3  to  0.6  m 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

>  0.6  m 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

Background  Alarm  Rates 

BAR*^”:  1  BAR*"”: 

Groups 

Found 

- 

- 

Identified 

- 

- 

Coverage 

-- 

- 

^In  cells  with  offset  data  entries,  the  numbers  to  the  left  are  the  result  and  the  two  numbers  to  the 
right  are  an  upper  and  lower  90-pereent  confidence  interval  for  an  assumed  binomial 
distribution. 

’’All  depths  are  measured  to  the  center  of  the  objeet. 
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TABLE  6f.  MOGUL  TEST  AREA  RESULTS  (not  covered) 


1  Response  Stage 

Discrimination  Stage  I 

Munitions'* 

Scores 

Small 

Medium 

Large 

Small 

Medium 

Large 

- 

- 

- 

- 

- 

- 

-- 

- 

1  By  Depth^  \ 

0  to  4D 

- 

- 

- 

- 

- 

- 

- 

- 

4D  to  8D 

- 

- 

- 

- 

- 

- 

- 

- 

8D  to  12D 

- 

- 

- 

- 

- 

- 

- 

- 

>  12D 

- 

- 

- 

- 

- 

- 

- 

- 

Clutter 

Scores 

Pfp 

1  By  Mass  \ 

By  Depth^ 

All 

Mass 

■m 

>  10  kg 

All 

Mass 

■m 

mm 

<  10kg 

All  Depth 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

0  to  0.15  m 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

0.15to0.3m 

- 

— 

- 

- 

— 

- 

- 

- 

— 

— 

0.3  to  0.6  m 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

>  0.6  m 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

Background  Alarm  Rates 

BAR™*:  1  BAR‘“*“: 

Groups 

Found 

- 

- 

Identifled 

- 

- 

Coverage 

- 

- 

^In  cells  with  offset  data  entries,  the  numbers  to  the  left  are  the  result  and  the  two  numbers  to  the 
right  are  an  upper  and  lower  90-pereent  eonfidenee  interval  for  an  assumed  binomial 
distribution. 

’’All  depths  are  measured  to  the  eenter  of  the  objeet. 


4,3  EFFICIENCY,  REJECTION  RATES,  AND  TYPE  CLASSIFICATION 

Effieieney  and  rejeetion  rates  are  ealeulated  to  quantify  the  diserimination  ability  at 
speeifie  points  of  interest  on  the  ROC  eurve:  (1)  at  the  point  where  no  deerease  in  is  suffered 
(i.e.,  the  effieieney  is  by  definition  equal  to  one)  and  (2)  at  the  operator  seleeted  threshold. 
These  values  are  presented  in  Tables  7a  through  7d. 


TABLE  7a,  BLIND  GRID  EFFICIENCY  AND 
REJECTION  RATES 


Efficiency  (E) 

False  Positive 
Rejection  Rate 

Background  Alarm 
Rejection  Rate 

At  Operating  Point 

0.96 

0.75 

0.82 

With  No  Loss  of  Pd 

1.00 

0.00 

0.00 
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TABLE  7b.  OPEN  FIELD  (DIRECT)  EFFICIENCY 
AND  REJECTION  RATES 


Efficiency  (E) 

False  Positive 
Rejection  Rate 

Background  Alarm 
Rejection  Rate 

At  Operating  Point 

0.97 

0.57 

0.68 

With  No  Loss  of  Pa 

1.00 

0.00 

0.00 

TABLE  7c.  OPEN  FIELD  (INDIRECT)  EFFICIENCY 
AND  REJECTION  RATES 


Efficiency  (E) 

False  Positive 
Rejection  Rate 

Background  Alarm 
Rejection  Rate 

At  Operating  Point 

0.94 

0.90 

0.92 

With  No  Loss  of  Pa 

1.00 

0.00 

0.00 

Note;  This  Results  Table  was  seored  using  a  subset  of  the  full  Open  Field  (Indireet  Fire) 
Ground  Truth,  therefore  these  results  cannot  be  directly  compared  to  any  other  APG 
Open  Field  (Indirect  Fire)  Results  reported. 


TABLE  7d.  OPEN  FIELD  (LEGACY)  EFFICIENCY  AND 
REJECTION  RATES 
(not  covered) 


Efficiency  (E) 

False  Positive 
Rejection  Rate 

Background  Alarm 
Rejion  Rate 

At  Operating  Point 

- 

- 

- 

With  No  Loss  of  Pa 

- 

- 

- 

TABLE  7e.  WOODED  EFFICIENCY  AND  REJECTION  RATES 

(not  covered) 


Efficiency  (E) 

False  Positive 
Rejection  Rate 

Background  Alarm 
Rejection  Rate 

At  Operating  Point 

- 

- 

- 

With  No  Loss  of  Pa 

- 

- 

- 
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TABLE  7f.  MOGUL  EFFICIENCY  AND  REJECTION  RATES 

(not  covered) 


Efficiency  (E) 

False  Positive 
Rejection  Rate 

Background  Alarm 
Rejection  Rate 

At  Operating  Point 

— 

— 

— 

With  No  Loss  of  Pd 

- 

- 

- 

At  the  demonstrator’s  reeommended  setting,  the  munitions  items  that  were  deteeted  and 
correctly  discriminated  were  further  scored  on  whether  their  correct  type  could  be  identified 
(table  8a  through  8f).  Correct  type  examples  include  20-mm  projectile,  105-mm  HEAT 
projectile,  and  2. 75 -inch  Rocket.  A  list  of  the  standard  type  declaration  required  for  each 
munitions  item  was  provided  to  demonstrators  prior  to  testing.  The  standard  types  for  the  three 
example  items  are  20-mmP,  105H,  and  2.75-inch. 


TABLE  8a,  BLIND  GRID  CORRECT  TYPE 
CLASSIFICATION  OF  TARGETS 
CORRECTLY  DISCRIMINATED 
AS  MUNITIONS 


Size 

Percentage  Correct,  % 

25mm 

40 

37mm 

100 

60mm 

87 

81mm 

73 

105mm 

13 

105  artillery 

100 

Overall 

69 

TABLE  8b.  OPEN  FIELD  DIRECT  FIRE 
CORRECT  TYPE  CLASSIFICATION 
OF  TARGETS  CORRECTLY 
DISCRIMINATED  AS 
MUNITIONS 


Size 

Percentage  Correct,  % 

25mm 

64 

37mm 

70 

105mm 

92 

Overall 

75 
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TABLE  8c.  OPEN  FIELD  INDIRECT  FIRE 
CORRECT  TYPE  CLASSIFICATION 
OF  TARGETS  CORRECTLY 
DISCRIMINATED  AS 
MUNITIONS 


Size 

Percentage  Correct,  % 

60  mm 

78 

81  mm 

69 

105  mm 

92 

Overall 

80 

Note;  This  Results  Table  was  seored  using  a  subset  of  the  full  Open  Field  (Indireet  Fire) 
Ground  Truth,  therefore  these  results  eannot  be  direetly  eompared  to  any  other  APG 
Open  Field  (Indireet  Fire)  Results  reported. 


TABLE  8d.  OPEN  FIELD  LEGACY  CORRECT 
TYPE  CLASSIFICATION  OF  TARGETS 
CORRECTLY  DISCRIMINATED 
AS  MUNITIONS  (not  covered) 


Size 

Percentage  Correct 

Small 

- 

Medium 

- 

Large 

- 

Overall 

- 

TABLE  8e.  WOODED  CORRECT  TYPE 
CLASSIFICATION  OF  TARGETS 
CORRECTLY  DISCRIMINATED 
AS  MUNITIONS  (not  covered) 


Size 

Percentage  Correct 

Small 

- 

Medium 

- 

Large 

- 

Overall 

- 
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TABLE  8f.  MOGUL  CORRECT  TYPE 
CLASSIFICATION  OF  TARGETS 
CORRECTLY  DISCRIMINATED 
AS  MUNITIONS  (not  covered) 


Size 

Percentage  Correct 

Small 

- 

Medium 

- 

Large 

- 

Overall 

- 

4,4  LOCATION  ACCURACY 

The  mean  loeation  error  and  standard  deviations  appear  in  Tables  9a  through  9f.  These 
calculations  are  based  on  average  missed  distance  for  munitions  correctly  identified  during  the 
response  stage.  Depths  are  measured  from  the  center  of  the  munitions  to  the  surface.  For  the 
blind  grid,  only  depth  errors  are  calculated  because  (X,  Y)  positions  are  known  to  be  the  centers 
of  the  grid  square. 


TABLE  9a,  BLIND  GRID  MEAN  LOCATION  ERROR 
AND  STANDARD  DEVIATION 


Mean 

standard  Deviation 

Northing 

NA 

NA 

Easting 

NA 

NA 

Depth 

0.00 

0.09 

TABLE  9b,  OPEN  FIELD  DIRECT  FIRE  MEAN 
LOCATION  ERROR  AND  STANDARD 
DEVIATION 


Mean 

Standard  Deviation 

Northing 

0.00 

0.09 

Easting 

-0.05 

0.10 

Depth 

0.00 

0.09 
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TABLE  9c.  OPEN  FIELD  INDIRECT  FIRE  MEAN  LOCATION 
ERROR  AND  STANDARD  DEVIATION 


Mean 

Standard  Deviation 

Northing 

0.01 

0.10 

Easting 

-0.05 

0.10 

Depth 

0.00 

0.09 

TABLE  9d.  OPEN  FIELD  LEGACY  MEAN  LOCATION 
ERROR  AND  STANDARD  DEVIATION  (not  covered) 


Mean 

Standard  Deviation 

Northing 

- 

- 

Easting 

- 

- 

Depth 

- 

- 

TABLE  9e.  WOODED  MEAN  LOCATION  ERROR 
AND  STANDARD  DEVIATION  (not  covered) 


Mean 

Standard  Deviation 

Northing 

- 

- 

Easting 

- 

- 

Depth 

- 

- 

TABLE  9f.  MOGUL  MEAN  LOCATION  ERROR 
AND  STANDARD  DEVIATION  (not  covered) 


Mean 

Standard  Deviation 

Northing 

- 

- 

Easting 

- 

- 

Depth 

- 

- 
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SECTIONS.  APPENDIXES 


APPENDIX  A,  TERMS  AND  DEFINITIONS 

GENERAL  DEEINITIONS 

Anomaly:  Location  of  a  system  response  deemed  to  warrant  further  investigation  by  the 
demonstrator  for  consideration  as  an  emplaced  munitions  item. 

Detection:  An  anomaly  location  that  is  within  Rhaio  of  an  emplaced  munitions  item. 

Military  Munitions  (MM):  Specific  categories  of  MM  that  may  pose  unique  explosive  safety 
risks,  including  UXO  as  defined  in  10  USC  101(e)(5),  DMM  as  defined  in  10  USC  2710(e)(2) 
and/or  munitions  constituents  (e.g.,  TNT,  RDX)  as  defined  in  10  USC  2710(e)(3)  that  are  present 
in  high  enough  concentrations  to  pose  an  explosive  hazard. 

Emplaced  Munitions:  A  munitions  item  buried  by  the  government  at  a  specified  location  in  the 
test  site. 

Emplaced  Clutter:  A  clutter  item  (i.e.,  nonmunitions  item)  buried  by  the  government  at  a 
specified  loeation  in  the  test  site. 

Rhaio:  A  predetermined  radius  about  an  emplaeed  item  (clutter  or  munitions)  within  which  an 
anomaly  identified  by  the  demonstrator  as  being  of  interest  is  eonsidered  to  be  a  deteetion  of  that 
item.  Eor  the  purpose  of  this  program,  a  circular  halo  0.5  meters  in  radius  is  placed  around  the 
eenter  of  the  object  for  all  clutter  and  munitions  items. 

Small  Munitions:  Caliber  of  munitions  less  than  or  equal  to  40  mm  (includes  20-mm  projeetile, 
25-mm  projectile,  37-mm  projectile,  40-mm  projectile,  submunitions  BLU-26,  BLU-63,  and 
M42). 

Medium  Munitions:  Caliber  of  munitions  greater  than  40  mm  and  less  than  or  equal  to  81  mm 
(includes  57-mm  projeetile,  60-mm  mortar,  2.75-inoh  roeket,  and  81 -mm  mortar). 

Large  Munitions:  Caliber  of  munitions  greater  than  81  mm  (ineludes  105 -mm  HEAT,  105 -mm 
projectile,  and  155-mm  projectile). 

Group:  Two  or  more  adjaeent  GT  items  with  overlapping  halos. 

GT :  Ground  truth 

Response  Stage  Noise  Level:  The  level  that  represents  the  signal  level  below  whieh  anomalies 
are  not  considered  detectable.  Demonstrators  are  required  to  provide  the  recommended  noise 
level  for  the  blind  grid  test  area. 


A-1 


Discrimination  Stage  Threshold;  The  demonstrator-seleeted  threshold  level  that  is  expeeted  to 
provide  optimum  performanee  of  the  system  by  retaining  all  deteetable  munitions  and  rejeeting 
the  maximum  amount  of  elutter.  This  level  defines  the  subset  of  anomalies  the  demonstrator 
would  reeommend  digging  based  on  diserimination. 

Binomially  Distributed  Random  Variable;  A  random  variable  of  the  type  whieh  has  only  two 
possible  outeomes,  say  sueeess  and  failure,  is  repeated  for  n  independent  trials  with  the 
probability  p  of  sueeess  and  the  probability  1-p  of  failure  being  the  same  for  eaeh  trial.  The 
number  of  sueeesses  x  observed  in  the  n  trials  is  an  estimate  of  p  and  is  eonsidered  to  be  a 
binomially  distributed  random  variable. 

RESPONSE  AND  DISCRIMINATION  STAGE  DATA 

The  seoring  of  the  demonstrator’s  performanee  is  eondueted  in  two  stages;  response  stage 
and  diserimination  stage.  Eor  both  stages,  the  probability  of  deteetion  (Pa)  and  the  false  alarms 
are  reported  as  reeeiver-operating  eharaeteristie  (ROC)  eurves.  Ealse  alarms  are  divided  into 
those  anomalies  that  eorrespond  to  emplaeed  elutter  items,  measuring  the  probability  of  elutter 
deteetion  (Pcd)  or  probability  of  false  positive  (Pfp).  Those  that  do  not  eorrespond  to  any  known 
item  are  termed  baekground  alarms. 

The  response  stage  is  a  measure  of  whether  the  sensor  ean  deteet  an  objeet  of  interest.  Eor 
a  ehannel  instrument,  this  value  should  be  elosely  related  to  the  amplitude  of  the  signal.  The 
demonstrator  must  report  the  response  level  (threshold)  below  whieh  target  responses  are 
deemed  insuffieient  to  warrant  further  investigation.  At  this  stage,  minimal  proeessing  may  be 
done.  This  ineludes  filtering  long-  and  short-seale  variations,  bias  removal,  and  sealing.  This 
proeessing  should  be  detailed  in  the  data  submission. 

Eor  a  multiehannel  instrument,  the  demonstrator  must  eonstruet  a  quantity  analogous  to 
amplitude.  The  demonstrator  should  eonsider  what  eombination  of  ehannels  provides  the  best 
test  for  deteeting  any  objeet  that  the  sensor  ean  deteet.  The  average  amplitude  aeross  a  set  of 
ehannels  is  an  example  of  an  aeeeptable  response  stage  quantity.  Other  methods  may  be  more 
appropriate  for  a  given  sensor.  Again,  minimal  proeessing  ean  be  done,  and  the  demonstrator 
should  explain  how  this  quantity  was  eonstrueted  in  their  data  submission. 

The  diserimination  stage  evaluates  the  demonstrator’s  ability  to  eorreetly  identify 
munitions  as  sueh,  and  to  rejeet  elutter.  Eor  the  same  loeations  as  in  the  response  stage  anomaly 
list,  the  diserimination  stage  list  eontains  the  output  of  the  algorithms  applied  in  the 
diserimination-stage  proeessing.  This  list  is  prioritized  based  on  the  demonstrator’s 
determination  that  an  anomaly  loeation  is  likely  to  eontain  munitions.  Thus,  higher  output  values 
are  indieative  of  higher  eonfidenee  that  a  munitions  item  is  present  at  the  speeified  loeation.  Eor 
eleetronie  signal  proeessing,  priority  ranking  is  based  on  algorithm  output.  Eor  other  systems, 
priority  ranking  is  based  on  human  judgment.  The  demonstrator  also  seleets  the  threshold  that 
the  demonstrator  believes  will  provide  optimum  system  performanee,  (i.e.,  that  retains  all  the 
deteeted  munitions  and  rejeets  the  maximum  amount  of  elutter). 

Note;  The  two  lists  provided  by  the  demonstrator  eontain  identieal  numbers  of  potential  target 
loeations.  They  differ  only  in  the  priority  ranking  of  the  deelarations. 
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GROUP  SCORING  FACTORS 


Based  on  configuration  of  the  GT  at  the  standardized  sites  and  the  defined  scoring 
methodology,  there  exists  munitions  groups  defined  as  having  overlapping  halos.  In  these  cases, 
the  following  scoring  logic  is  implemented  (fig.  A-1  through  A-9); 

a.  Overall  site  scores  (i.e.,  Pd)  will  consider  only  isolated  munitions  and  clutter  items. 

b.  GT  items  that  have  overlapping  halos  (both  munitions  and  clutter)  will  form  a  group 
and  groups  may  form  chains. 

c.  Groups  will  have  a  complex  halos  composed  of  all  the  composite  halos  of  all  its  GT 

items. 

d.  Groups  will  have  three  scoring  factors:  groups  found  groups  identified  and  group 
coverage.  Scores  will  be  based  on  1 : 1  matches  of  anomalies  and  GT. 

(1)  Groups  Found  (Found):  the  number  of  groups  that  have  one  or  more  GT  items 
matched  divided  by  the  total  number  of  groups.  Demonstrators  will  be  credited  with  detecting  a 
group  if  any  item  within  the  group  is  matched  to  an  anomaly  in  their  list. 

(2)  Groups  Identified  (ID):  the  number  of  groups  that  have  two  or  more  GT  items 
matched  divided  by  the  total  number  of  groups.  Demonstrators  will  be  credited  with  identifying 
that  a  group  is  present  if  multiple  items  within  the  composite  halo  are  matched  to  anomalies  in 
their  list. 

(3)  Group  Coverage  (Coverage):  the  number  of  GT  items  matched  within  groups  divided 
by  the  total  number  of  GT  items  within  groups.  This  metric  measures  the  demonstrator  accuracy 
in  determining  the  number  of  anomalies  within  a  group.  If  five  items  are  present  and  only  two 
anomalies  are  matched,  the  demonstrator  will  score  0.4.  If  all  five  are  matched  the  demonstrator 
will  score  1.0. 

e.  Location  error  will  not  be  reported  for  groups. 
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f.  Demonstrators  will  not  be  asked  to  call  out  groups  in  their  scoring  submissions.  If 
multiple  anomalies  are  indicated  in  a  small  area,  the  demonstrator  will  report  all  individual 
anomalies. 

g.  Excess  alarms  within  a  halo  will  be  disregarded. 


A-1.  Example  of  detected  item. 
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A-3.  Example  of  group  identified  (ID). 


A-4.  Example  of  excess  alarms  disregarded. 
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A-6.  Example  of  group  (1/4  =  0.25). 
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A-7.  Example  of  group  (2/4  =  0.5). 
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RESPONSE  STAGE  DEEINITIONS 

Response  Stage  Probability  of  Detection  (P/®’^):  =  (No.  of  response-stage  detections)/ 

(No.  of  emplaced  munitions  in  the  test  site). 

Response  Stage  Clutter  Detection  (cd’^®®):  An  anomaly  location  that  is  within  Rhaio  of  an 
emplaced  clutter  item. 

Response  Stage  Probability  of  Clutter  Detection  (Pc/'^'*):  =  (No.  of  response-stage  clutter 

detections)/(No.  of  emplaced  clutter  items). 

Response  Stage  Background  Alarm  (ba’^^^®):  An  anomaly  in  a  blind  grid  cell  that  contains  neither 
emplaced  munitions  nor  an  emplaced  clutter  item.  An  anomaly  location  in  the  open  field  or 
scenarios  that  is  outside  Rhaio  of  any  emplaced  munitions  or  emplaced  clutter  item. 

Response  Stage  Probability  of  Background  Alarm  (Pba’^^’^):  Blind  grid  only:  Pba'^'^’^  =  (No.  of 
response-stage  background  alarms)/(No.  of  empty  grid  locations). 

Response  Stage  Background  Alarm  Rate  (BAR’^'^'*):  Open  field  any  challenge  area  (including  the 
direct  and  indirect  firing  sub  areas)  only:  BAR’^^^'*  =  (No.  of  response-stage  background 
alarms)/(arbitrary  constant). 

Note  that  the  quantities  P/®^  Pc/^^  Pba‘^^^  and  BAR’^®’^  are  functions  of  the  threshold 
applied  to  the  response-stage  signal  strength.  These  quantities  can  therefore  be  written  as 
PdreS(tres)^  Pe/®\f®*),  Pb/“(f''),  and  B AR^^^f®^) . 
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DISCRIMINATION  STAGE  DEFINITIONS 


Discrimination;  The  application  of  a  signal  processing  algorithm  or  human  judgment  to  sensor 
data  to  diseriminate  munitions  from  clutter.  Diserimination  should  identify  anomalies  that  the 
demonstrator  has  high  confidence  correspond  to  munitions,  as  well  as  those  that  the  demonstrator 
has  high  eonfidenee  eorrespond  to  nonmunitions  or  background  returns.  The  former  should  be 
ranked  with  highest  priority  and  the  latter  with  lowest. 

Diserimination  Stage  Probability  of  Detection  (P/'“):  =  (No.  of  discrimination-stage 

detections)/(No.  of  emplaced  munitions  in  the  test  site). 

Diserimination  Stage  False  Positive  (fp‘^‘®°):  An  anomaly  loeation  that  is  within  Rhaio  of  an 
emplaced  elutter  item. 

Diserimination  Stage  Probability  of  False  Positive  (Pfp‘*''^‘^):  =  (No.  of  diserimination  stage 

false  positives)/(No.  of  emplaced  clutter  items). 

Diserimination  Stage  Background  Alarm  (ba‘^‘®‘^):  An  anomaly  in  a  blind  grid  cell  that  eontains 
neither  emplaced  munitions  nor  an  emplaeed  clutter  item.  An  anomaly  loeation  in  the  open  field 
or  scenarios  that  is  outside  Rhaio  of  any  emplaced  munitions  or  emplaced  elutter  item. 

Diserimination  Stage  Probability  of  Background  Alarm  (Pba^'’^^^):  Pba^'*°  =  (No.  of  discrimination- 
stage  background  alarms)/(No.  of  empty  grid  locations). 

Discrimination  Stage  Background  Alarm  Rate  (BAR‘^‘®‘^):  BAR‘^‘®‘^  =  (No.  of  diserimination-stage 
baekground  alarms)/(arbitrary  constant). 

Note  that  the  quantities  Pfp‘*'*“,  Pba‘^‘*‘^,  and  BAR‘*'®‘^  are  functions  of  t'*'®®,  the  threshold 
applied  to  the  diserimination-stage  signal  strength.  These  quantities  ean  therefore  be  written  as 

PdSisC(td:sC)^  p^d:sc(tdisc)^  Pba‘^‘'‘^(t‘*'^‘=),  and  B AR‘^‘*‘=(t‘^‘^‘=) . 


A-9 


RECEIVER-OPERATING  CHARACTERISTIC  (ROC)  CURVES 


ROC  curves  at  both  the  response  and  diserimination  stages  ean  be  constructed  based  on  the 
above  definitions.  The  ROC  curves  plot  the  relationship  between  Pd  versus  Pcd  or  Pfj,  and  Pd 
versus  BAR  or  Pba  as  the  threshold  applied  to  the  signal  strength  is  varied  from  its  minimum 
(tmin)  to  its  maximum  (tmax)  value. '  Pd  versus  Pfp  and  Pd  versus  BAR  being  combined  into  ROC 
curves  are  shown  in  Eigure  A- 10.  Note  that  the  “res”  and  “disc”  superscripts  have  been 
suppressed  from  all  the  variables  for  clarity. 


Eigure  A- 10.  ROC  curves  for  open  field  testing.  Eaeh  curve  applies  to  both  the  response  and 
diserimination  stages. 

METRICS  TO  CHARACTERIZE  THE  DISCRIMINATION  STAGE 

The  demonstrator  is  also  scored  on  efficieney  and  rejection  ratio,  which  measure  the 
effeetiveness  of  the  discrimination  stage  proeessing.  The  goal  of  discrimination  is  to  retain  the 
greatest  number  of  munitions  detections  from  the  anomaly  list  while  rejecting  the  maximum 
number  of  anomalies  arising  from  nonmunitions  items.  The  efficiency  measures  the  fraction  of 
detected  munitions  retained  by  the  diserimination,  while  the  rejection  ratio  measures  the  fraction 
of  false  alarms  rejected.  Both  measures  are  defined  relative  to  the  entire  response  list,  i.e.,  the 
maximum  munitions  detectable  by  the  sensor  and  its  aecompanying  elutter  detection  rate/false 
positive  rate  or  baekground  alarm  rate. 

'Strictly  speaking,  ROC  curves  plot  the  Pd  versus  Pba  over  a  predetermined  and  fixed  number  of 
deteetion  opportunities  (some  of  the  opportunities  are  located  over  munitions  and  others  are 
loeated  over  elutter  or  blank  spots).  In  an  open  field  scenario,  each  system  suppresses  its  signal 
strength  reports  until  some  bare -minimum  signal  response  is  received  by  the  system. 
Consequently,  the  open  field  ROC  curves  do  not  have  information  from  low  signal-output 
loeations,  and,  furthermore,  different  eontraetors  report  their  signals  over  a  different  set  of 
locations  on  the  ground.  These  ROC  curves  are  thus  not  true  to  the  strict  definition  of  ROC 
curves  as  defined  in  textbooks  on  detection  theory.  Note,  however,  that  the  ROC  eurves 
obtained  in  the  blind  grid  test  sites  are  true  ROC  eurves. 
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Efficiency  (E):  E  =  Measures  (at  a  threshold  of  interest)  the  degree 

to  whieh  the  maximum  theoretieal  detection  performance  of  the  sensor  system  (as  determined  by 
the  response  stage  tmin)  is  preserved  after  applieation  of  diserimination  teehniques.  Efficieney  is 
a  number  between  0  and  1 .  An  efficiency  of  1  implies  that  all  of  the  munitions  initially  detected 
in  the  response  stage  were  retained  at  the  specified  threshold  in  the  diserimination  stage, 

Ealse  Positive  Rejection  Rate  (Rfp):  Rfj,  =  1  -  [Pfp‘*'*‘^(t‘^‘^‘^)/Pc/^*(tmin'^^'^)]:  Measures  (at  a 
threshold  of  interest)  the  degree  to  which  the  sensor  system's  false  positive  performance  is 
improved  over  the  maximum  false  positive  performanee  (as  determined  by  the  response  stage 
tmin).  The  rejection  rate  is  a  number  between  0  and  E  A  rejeetion  rate  of  1  implies  that  all 
emplaced  clutter  initially  detected  in  the  response  stage  were  correctly  rejected  at  the  specified 
threshold  in  the  diserimination  stage. 

Baekground  Alarm  Rejection  Rate  (Rba): 

Blind  grid:  Rba  =  1  -  [Pba"'^‘^(t""^)/Pba''^\t^,n''‘^^)]. 

Open  field:  Rba  =  1  -  [BAR‘*'"‘=(t‘^'"‘^)/BAR^®*(Wn''®")]). 

Measures  the  degree  to  whieh  the  diserimination  stage  correctly  rejects  background  alarms 
initially  detected  in  the  response  stage.  The  rejection  rate  is  a  number  between  0  and  1.  A 
rejection  rate  of  1  implies  that  all  baekground  alarms  initially  deteeted  in  the  response  stage  were 
rejected  at  the  speeified  threshold  in  the  discrimination  stage. 

CHI-SQUARE  COMPARISON 

The  Chi-square  test  for  differences  in  probabilities  (or  2  by  2  contingency  table)  is  used  to 
analyze  two  samples  drawn  from  two  different  populations  to  see  if  both  populations  have  the 
same  or  different  proportions  of  elements  in  a  eertain  category.  More  specifieally,  two  random 
samples  are  drawn,  one  from  eaeh  population,  to  test  the  null  hypothesis  that  the  probability  of 
event  A  (some  specified  event)  is  the  same  for  both  populations. 

The  test  statistic  of  the  2  by  2  eontingency  table  is  the  Chi-square  distribution  with  one 
degree  of  freedom.  When  an  assoeiation  between  a  more  ehallenging  terrain  feature  and 
relatively  degraded  performance  is  sought,  a  one-sided  test  is  performed.  A  two-sided  2  by 
2  eontingency  table  is  used  in  the  Standardized  UXO  Teehnology  Demonstration  Site  Program  to 
eompare  performanee  between  any  two  areas  or  subareas  when  the  direetion  of  degradation 
cannot  be  predetermined. 

Eor  a  one-sided  test,  a  significance  level  of  0.05  is  used  to  set  the  eritieal  deeision  limit.  It 
is  a  eritieal  decision  limit  because  if  the  test  statistie  ealculated  from  the  data  exceeds  this  value, 
then  the  lower  proportion  tested  will  be  eonsidered  signifieantly  less  than  the  greater  one 
(degraded).  If  the  test  statistic  calculated  from  the  data  is  less  than  this  value,  then  no 
degradation  can  be  said  to  exist  beeause  of  the  terrain  feature  introdueed. 
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For  a  two-sided  test,  a  signifieanee  level  of  0.10  is  used  to  allow  .05  on  either  side  of  the 
deeision.  It  is  a  eritieal  deeision  limit  beeause  if  the  test  statistie  ealculated  from  the  data 
exeeeds  this  value,  then  the  two  proportions  tested  will  be  eonsidered  signifieantly  different.  If 
the  test  statistic  calculated  from  the  data  is  less  than  this  value,  then  the  two  proportions  tested 
will  be  considered  not  significantly  different. 

An  exception  must  be  applied  when  either  a  0  or  100  percent  success  rate  occurs  in  the 
sample  data.  The  Chi-square  test  cannot  be  used  in  these  instances.  Instead,  Fischer’s  test  is 
used,  and  the  critical  decision  limit  for  one-sided  tests  is  the  chosen  significance  level,  which  in 
this  case  is  0.05.  With  Fischer’s  test,  if  the  test  statistic  is  less  than  the  critical  value,  then  the 
proportions  are  considered  to  be  significantly  different. 

An  example  follows  that  illustrates  Standardized  UXO  Technology  Demonstration  Site 
blind  grid  results  compared  to  those  from  the  open  field  legacy.  It  should  be  noted  that  a 
significant  result  does  not  prove  a  cause-and-effect  relationship  exists  between  the  two 
populations  of  interest;  however,  it  does  serve  as  a  tool  to  indicate  that  one  data  set  has 
experienced  a  degradation  or  change  in  system  performance  at  a  large  enough  level  than  can  be 
accounted  for  merely  by  chance  or  random  variation.  Note  also  that  a  result  that  is  not 
significant  indicates  that  there  is  not  enough  evidence  to  declare  that  anything  more  than  chance 
or  random  variation  within  the  same  population  is  at  work  between  the  two  data  sets  being 
compared. 

Demonstrator  X  achieves  the  following  overall  results  after  surveying  the  blind  grid  and 
open  field  (legacy)  using  the  same  system  (results  indicate  the  number  of  munitions  detected 
divided  by  the  number  of  munitions  emplaced): 

Blind  grid  Open  field 

P/"M00/100  =  1.0  8/10  =  .80 


P/®*:  BLIND  GRID  versus  OPEN  FIELD  (legacy).  Using  the  example  data  above  to 
compare  probabilities  of  detection  in  the  response  stage,  all  100  munitions  out  of  100  emplaced 
munitions  items  were  detected  in  the  blind  grid  while  8  munitions  out  of  10  emplaced  were 
detected  in  the  open  field.  Eischer’s  test  must  be  used  since  a  100  percent  success  rate  occurs  in 
the  data.  Eischer’s  test  uses  the  four  input  values  to  calculate  a  test  statistic  of  0.0075  that  is 
compared  against  the  critical  value  of  0.05.  Since  the  test  statistic  is  less  than  the  critical  value, 
the  smaller  response  stage  detection  rate  (0.80)  is  considered  to  be  significantly  less  at  the 
0.05  level  of  significance.  While  a  significant  result  does  not  prove  a  cause-and-effect 
relationship  exists  between  the  change  in  survey  area  and  degradation  in  performance,  it  does 
indicate  that  the  detection  ability  of  demonstrator  X’s  system  seems  to  have  been  degraded  in  the 
open  field  relative  to  results  from  the  blind  grid  using  the  same  system.  This  is  an  example  of  a 
one-sided  Chi-squared  test. 
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APPENDIX  B.  DAILY  WEATHER  LOGS 


Date,  10 

Time,  “EST 

Avg.  Temp,  °F 

Total  Precip.,  in. 

2  Mar 

0700 

33.3 

0.00 

0800 

34.2 

0.00 

0900 

37.2 

0.00 

1000 

40.8 

0.00 

1100 

42.3 

0.00 

1200 

43.0 

0.00 

1300 

43.3 

0.00 

1400 

43.7 

0.00 

1500 

43.5 

0.00 

1600 

42.6 

0.00 

1700 

41.0 

0.00 

3  Mar 

0700 

37.8 

0.00 

0800 

38.1 

0.00 

0900 

38.5 

0.00 

1000 

39.0 

0.00 

1100 

39.6 

0.00 

1200 

40.1 

0.00 

1300 

41.4 

0.00 

1400 

41.7 

0.00 

1500 

42.6 

0.00 

1600 

42.8 

0.00 

1700 

43.0 

0.00 

4  Mar 

0700 

35.6 

0.00 

0800 

36.1 

0.00 

0900 

38.1 

0.00 

1000 

40.5 

0.00 

1100 

41.9 

0.00 

1200 

42.1 

0.00 

1300 

43.2 

0.00 

1400 

44.8 

0.00 

1500 

45.7 

0.00 

1600 

45.9 

0.00 

1700 

45.5 

0.00 

“Eastern  Standard  Time. 


B-1 


Date,  10 

Time,  "EST 

Avg.  Temp,  °F 

Total  Precip.,  in. 

5  Mar 

0700 

36.5 

0.00 

0800 

36.7 

0.00 

0900 

38.5 

0.00 

1000 

39.7 

0.00 

1100 

41.2 

0.00 

1200 

44.4 

0.00 

1300 

46.4 

0.00 

1400 

47.8 

0.00 

1500 

47.3 

0.00 

1600 

48.0 

0.00 

1700 

47.5 

0.00 

8  Mar 

0700 

35.6 

0.00 

0800 

39.0 

0.00 

0900 

44.2 

0.00 

1000 

49.8 

0.00 

1100 

52.5 

0.00 

1200 

54.5 

0.00 

1300 

56.7 

0.00 

1400 

58.3 

0.00 

1500 

59.7 

0.00 

1600 

60.1 

0.00 

1700 

59.9 

0.00 

9  Mar 

0700 

36.7 

0.00 

0800 

39.0 

0.00 

0900 

50.7 

0.00 

1000 

54.3 

0.00 

1100 

57.2 

0.00 

1200 

59.5 

0.00 

1300 

60.8 

0.00 

1400 

61.2 

0.00 

1500 

62.1 

0.00 

1600 

62.2 

0.00 

1700 

61.5 

0.00 

'Eastern  Standard  Time. 


APPENDIX  C.  SOIL  MOISTURE 


Date:  3  Mar  10 

Time:  NA,  1400 

Probe  Location 

Layer,  in. 

A,M,  Reading,  % 

P,M,  Reading,  % 

Wet  area 

0  to  6 

— 

— 

6  to  12 

— 

— 

12  to  24 

— 

— 

24  to  36 

— 

— 

36  to  48 

— 

— 

Wooded  area 

0  to  6 

— 

— 

6  to  12 

— 

— 

12  to  24 

— 

— 

24  to  36 

— 

— 

36  to  48 

— 

— 

Open  area 

0  to  6 

— 

— 

6  to  12 

— 

— 

12  to  24 

— 

— 

24  to  36 

— 

— 

36  to  48 

— 

— 

Calibration  lanes 

0  to  6 

— 

28.9 

6  to  12 

— 

31.7 

12  to  24 

— 

34.8 

24  to  36 

— 

39.7 

36  to  48 

— 

44.1 

Blind  grid/moguls 

0  to  6 

— 

— 

6  to  12 

— 

— 

12  to  24 

— 

— 

24  to  36 

— 

— 

36  to  48 

- 

- 

C-1 


C-2 


C-3 


C-4 


C-5 


C-6 


Status 

Stop 

Duration 

Operational  Status 

Track 

Date 

■mIHB 

AreaTested 

Time 

min. 

Operational  Status 

Comments 

Method 

Pattern 

Field  Conditions  | 

3/2/2010 

3 

CALIBRATION  LANES 

1445 

1600 

75 

INITIAL  SET-UP 

INITIAL  SET  UP 

GPS 

LINEAR 

RAINY 

MUDDY 

3/3/2010 

3 

CALIBRATION  LANES 

730 

1155 

265 

INITIAL  SET-UP 

INITIAL  SET  UP 

GPS 

LINEAR 

RAINY 

MUDDY 

3/3/2010 

3 

CALIBRATION  LANES 

1155 

1425 

150 

COLLECTING 

DATA 

COLLECT  DATA 

GPS 

LINEAR 

RAINY 

MUDDY 

3/3/2010 

3 

CALIBRATION  LANES 

1425 

1440 

15 

CALIBRATION 

CALIBRATE 

GPS 

LINEAR 

RAINY 

MUDDY 

3/3/2010 

3 

CALIBRATION  LANES 

1440 

1500 

20 

DOWNTIME  DUE 

TO  EQUIP 
MAINT/CHECK 

DOWNLOAD 

DATA 

GPS 

LINEAR 

RAINY 

MUDDY 

3/3/2010 

3 

CALIBRATION  LANES 

1500 

1515 

15 

BREAK/LUNCH 

BREAK/LUNCH 

GPS 

LINEAR 

RAINY 

MUDDY 

DOWNTIME  DUE 

3/3/2010 

3 

CALIBRATION  LANES 

1515 

1540 

25 

TO  EQUIP 
MAINT/CHECK 

DATA  CHECK 

GPS 

LINEAR 

RAINY 

MUDDY 

3/3/2010 

3 

CALIBRATION  LANES 

1540 

1600 

20 

DAILY  START, 

STOP 

EQUIPMENT 

BREAKDOWN 

GPS 

LINEAR 

RAINY 

MUDDY 

3/4/2010 

3 

BLIND  TEST  GRID 

740 

825 

45 

DAILY  START, 

STOP 

SETUP 

EQUIPMENT 

GPS 

LINEAR 

SUNNY 

MUDDY 

3/4/2010 

3 

BLIND  TEST  GRID 

825 

845 

20 

CALIBRATION 

CALIBRATE 

GPS 

LINEAR 

SUNNY 

MUDDY 

3/4/2010 

3 

BLIND  TEST  GRID 

845 

1235 

230 

COLLECTING 

DATA 

COLLECT  DATA 

GPS 

LINEAR 

SUNNY 

MUDDY 

3/4/2010 

3 

BLIND  TEST  GRID 

1235 

1300 

25 

CALIBRATION 

CALIBRATE 

GPS 

LINEAR 

SUNNY 

MUDDY 

3/4/2010 

3 

BLIND  TEST  GRID 

1300 

1335 

35 

BREAK/LUNCH 

BREAK/LUNCH 

GPS 

LINEAR 

SUNNY 

MUDDY 

3/4/2010 

3 

OPEN  FIELD 

1335 

1450 

75 

DAILY  START, 

STOP 

SETUP 

EQUIPMENT 

GPS 

LINEAR 

SUNNY 

MUDDY 

3/4/2010 

3 

OPEN  FIELD 

1450 

1615 

85 

COLLECTING 

DATA 

COLLECT  DATA, 
INDIRECT  FIRE 

GPS 

LINEAR 

SUNNY 

MUDDY 

3/4/2010 

3 

OPEN  FIELD 

1615 

1635 

20 

CALIBRATION 

CALIBRATE 

GPS 

LINEAR 

SUNNY 

MUDDY 

3/4/2010 

3 

OPEN  FIELD 

1635 

1700 

25 

DAILY  START, 

STOP 

EQUIPMENT 

BREAKDOWN 

GPS 

LINEAR 

SUNNY 

MUDDY 

3/5/2010 

3 

OPEN  FIELD 

735 

905 

90 

DAILY  START, 

STOP 

SETUP 

EQUIPMENT 

GPS 

LINEAR 

SUNNY 

MUDDY 

3/5/2010 

3 

OPEN  FIELD 

905 

920 

15 

CALIBRATION 

CALIBRATE 

GPS 

LINEAR 

SUNNY 

MUDDY 

3/5/2010 

3 

OPEN  FIELD 

920 

1020 

60 

DOWNTIME  DUE 

TO  EQUIPMENT 
FAILURE 

FALUTY  GPS 
CABLE,  SODERED 

GPS 

LINEAR 

SUNNY 

MUDDY 

3/5/2010 

3 

OPEN  FIELD 

1020 

1030 

10 

CALIBRATION 

CALIBRATE 

GPS 

LINEAR 

SUNNY 

MUDDY 

3/5/2010 

3 

OPEN  FIELD 

1030 

1200 

90 

COLLECTING 

DATA 

COLLECT  DATA, 
INDIRECT  FIRE 

GPS 

LINEAR 

SUNNY 

MUDDY 

APPENDIX  D.  DAILY  ACTIVITY  LOGS 


D-2 


Date 


3/5/2010 


3/5/2010 


3/5/2010 


3/5/2010 


3/5/2010 


AreaTested 


OPEN  FIELD 


OPEN  FIELD 


OPEN  FIELD 


OPEN  FIELD 


OPEN  FIELD 


Duration 

min. 


Operational  Status 


DAILY  START, 
STOP 


BREAK/LUNCH 


COLLECTING 

DATA 


CALIBRATION 


DAILY  START, 
STOP 


Operational  Status 
Comments 


SETUP 

EQUIPMENT 


BREAK/LUNCH 


COLLECT  DATA, 
INDIRECT  FIRE 


CALIBRATE 


EQUIPMENT 

BREAKDOWN 


Track 

Method  Pattern 


GPS  LINEAR 


GPS  LINEAR 


GPS  LINEAR 


GPS  LINEAR 


GPS  LINEAR 


Field  Conditions 


SUNNY  MUDDY 


SUNNY  MUDDY 


SUNNY  MUDDY 


SUNNY  MUDDY 


SUNNY  MUDDY 


3/8/2010 

3 

OPEN  FIELD 

735 

800 

25 

DAILY  START, 

STOP 

SETUP 

EQUIPMENT 

GPS 

LINEAR 

SUNNY 

MUDDY 

3/8/2010 

3 

OPEN  FIELD 

800 

830 

30 

CALIBRATION 

CALIBRATE 

GPS 

LINEAR 

SUNNY 

MUDDY 

3/8/2010 

3 

OPEN  FIELD 

830 

1055 

145 

COLLECTING 

DATA 

COLLECT  DATA, 
INDIRECT  FIRE 

GPS 

LINEAR 

SUNNY 

MUDDY 

3/8/2010 

3 

OPEN  FIELD 

1055 

1125 

30 

DOWNTIME  DUE 

TO  EQUIPMENT 
FAILURE 

FLAT  TIRE,  FIXED 

GPS 

LINEAR 

SUNNY 

MUDDY 

3/8/2010 

3 

OPEN  FIELD 

1125 

1530 

245 

COLLECTING 

DATA 

COLLECT  DATA, 
INDIRECT  FIRE 

GPS 

LINEAR 

SUNNY 

MUDDY 

3/8/2010 

3 

OPEN  FIELD 

1530 

1540 

10 

DOWNTIME  DUE 

TO  EQUIPMENT 
FAILURE 

FLAT  TIRE,  FIXED 

GPS 

LINEAR 

SUNNY 

MUDDY 

3/8/2010 

3 

OPEN  FIELD 

1540 

1625 

45 

COLLECTING 

DATA 

COLLECT  DATA, 
INDIRECT  FIRE 

GPS 

LINEAR 

SUNNY 

MUDDY 

3/8/2010 

3 

OPEN  FIELD 

1625 

1635 

10 

CALIBRATION 

CALIBRATE 

GPS 

LINEAR 

SUNNY 

MUDDY 

3/8/2010 

3 

OPEN  FIELD 

1635 

1700 

25 

DAILY  START, 

STOP 

EQUIPMENT 

BREAKDOWN 

GPS 

LINEAR 

SUNNY 

MUDDY 

3/9/2010 

3 

OPEN  FIELD 

735 

805 

30 

DAILY  START, 

STOP 

SETUP 

EQUIPMENT 

GPS 

LINEAR 

SUNNY 

MUDDY 

3/9/2010 

3 

OPEN  FIELD 

805 

835 

30 

CALIBRATION 

CALIBRATE 

GPS 

LINEAR 

SUNNY 

MUDDY 

3/9/2010 

3 

OPEN  FIELD 

835 

1115 

160 

COLLECTING 

DATA 

COLLECT  DATA, 
INDIRECT  FIRE 

GPS 

LINEAR 

SUNNY 

MUDDY 

3/9/2010 

3 

OPEN  FIELD 

1115 

1135 

20 

DAILY  START, 

STOP 

SETUP 

EQUIPMENT 

GPS 

LINEAR 

SUNNY 

MUDDY 

3/9/2010 

3 

OPEN  FIELD 

1135 

1155 

20 

BREAK/LUNCH 

BREAK/LUNCH 

GPS 

LINEAR 

SUNNY 

MUDDY 

3/9/2010 

3 

OPEN  FIELD 

1155 

1550 

235 

COLLECTING 

DATA 

COLLECT  DATA, 
DIRECT  FIRE 

GPS 

LINEAR 

SUNNY 

MUDDY 

3/9/2010 

3 

OPEN  FIELD 

1550 

1555 

5 

DOWNTIME  DUE 

TO  EQUIP 
MAINT/CHECK 

CHANGE 

BATTERIES 

GPS 

LINEAR 

SUNNY 

MUDDY 

3/9/2010 

3 

OPEN  FIELD 

1555 

1630 

35 

COLLECTING 

DATA 

COLLECT  DATA, 
DIRECT  FIRE 

GPS 

LINEAR 

SUNNY 

MUDDY 

Date 

No.  of 
People 

Area  Tested 

Status 

Start 

Time 

Status 

Stop 

Time 

Duration 

min. 

Operational  Status 

Operational  Status 
Comments 

Track 

Method 

Pattern 

Field  Conditions 

3/9/2010 

3 

OPEN  FIELD 

1630 

1640 

10 

CALIBRATION 

CALIBRATE 

GPS 

LINEAR 

SUNNY 

MUDDY 

3/9/2010 

3 

OPEN  FIELD 

1640 

1700 

20 

DAILY  START, 

STOP 

EQUIPMENT 

BREAKDOWN 

GPS 

LINEAR 

SUNNY 

MUDDY 

3/10/2010 

3 

OPEN  FIELD 

740 

800 

20 

DAILY  START, 

STOP 

SETUP 

EQUIPMENT 

GPS 

LINEAR 

SUNNY 

MUDDY 

3/10/2010 

3 

OPEN  FIELD 

800 

835 

35 

CALIBRATION 

CALIBRATE 

GPS 

LINEAR 

SUNNY 

MUDDY 

3/10/2010 

3 

OPEN  FIELD 

835 

1205 

210 

COLLECTING 

DATA 

COLLECT  DATA, 
DIRECT  FIRE 

GPS 

LINEAR 

SUNNY 

MUDDY 

3/10/2010 

3 

OPEN  FIELD 

1205 

1220 

15 

COLLECTING 

DATA 

COLLECT  DATA, 
INDIRECT  FIRE 

GPS 

LINEAR 

SUNNY 

MUDDY 

3/10/2010 

3 

OPEN  FIELD 

1220 

1235 

15 

CALIBRATION 

CALIBRATE 

GPS 

LINEAR 

SUNNY 

MUDDY 

3/10/2010 

3 

OPEN  FIELD 

1235 

1320 

45 

DOWNTIME  DUE 

TO  EQUIP 
MAINT/CHECK 

DOWNLOAD 

DATA 

GPS 

LINEAR 

SUNNY 

MUDDY 

3/10/2010 

3 

OPEN  FIELD 

1320 

1630 

190 

DEMOBILIZATION 

DEMOBILIZATION 

GPS 

LINEAR 

SUNNY 

MUDDY 

APPENDIX  E,  REFERENCES 


1 .  Standardized  UXO  Technology  Demonstration  Site  Handbook,  DTC  Project 
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2.  Aberdeen  Proving  Ground  Soil  Survey  Report,  October  1998. 

3.  Data  Summary,  UXO  Standardized  Test  Site:  APG  Soils  Description,  May  2002. 
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APPENDIX  F.  ABBREVIATIONS 


APG  =  Aberdeen  Proving  Ground 

ATC  =  U.S.  Army  Aberdeen  Test  Center 

ATSS  =  Aberdeen  Test  Support  Services 

BAR  =  background  alarm  rate 

DMM  =  discarded  military  munitions 

EMI  =  electromagnetic  interference 

EQT  =  Environmental  Quality  Technology 

ERDC  =  U.S.  Army  Corps  of  Engineers  Engineering  Research  and 
Development  Center 
EST  =  Eastern  Standard  Time 

ESTCP  =  Environmental  Security  Technology  Certification  Program 
GPS  =  Global  Positioning  System 

GT  =  ground  truth 

HDSD  =  Homeland  Defense  and  Sustainment  Division 

HEAT  =  high-explosive  antitank 

JPG  =  Jefferson  Proving  Ground 

MM  =  military  munitions 

NS  =  nonstandard  munition 

POC  =  point  of  contact 

QA  =  quality  assurance 

QC  =  quality  control 

ROC  =  receiver-operating  characteristic 

S  =  standard  munition 

SAIC  =  Science  Applications  International  Corporation 

SCEMP  =  Simplified  Combined  EMI  Magnetometer  Prototype 

SERDP  =  Strategic  Environmental  Research  and  Development  Program 

USAEC  =  U.S.  Army  Environmental  Command 

UXO  =  unexploded  ordnance 

YPG  =  U.S.  Army  Yuma  Proving  Ground 
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APPENDIX  G.  DISTRIBUTION  LIST 


DTC  Project  No.  8-CO-160-UXO-021 


Addressee 


Commander 

U.S.  Army  Aberdeen  Test  Center 

ATTN;  TEDT-AT-SLE  (Mr.  J.  Stephen  McClung) 

400  Colleran  Road 

Aberdeen  Proving  Ground,  MD  21005-5059 

Program  Manager 
SERDP/ESTCP 
Munitions  Management 
ATTN:  Mr.  Herb  Nelson 
901  North  Stuart  Street,  Suite  303 
Arlington,  VA  22203 

U.S.  Geological  Survey 
ATTN:  Mr.  TedAsch 
Denver  Eederal  Center 
Bldg.  20,  MS-964 
Denver,  CO  80225-0046 

Defense  Technical  Information  Center 
8725  John  J.  Kingman  Road,  Suite  0944 
Port  Belvoir,  VA  22060-6218 

U.S.  Army  Environmental  Command 
ATTN:  IMAE-IT  (Dr.  Robert  Kirgan,  K-16) 

1711  1H35,  Suite  110 
San  Antonio,  TX  78233 


Secondary  distribution  is  controlled  by  Program  Manager,  SERDP/ESTCP,  Munitions 
Management,  ATTN;  Mr.  Herb  Nelson. 


No.  of 
Copies 


1 


1 


1 


PDP 


PDP 
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